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Abstract 

Path integration is a navigation strategy widely observed in nature where an 
animal maintains a running estimate, called the home vector, of its location during 
an excursion. Evidence suggests it is both ancient and ubiquitous in nature, and has 
been studied for over a century. In that time, canonical and neural network models 
have flourished, based on a wide range of assumptions, justifications and supporting 
data. Despite the importance of the phenomenon, consensus and unifying principles 
appear lacking. A fundamental issue is the neural representation of space needed for 
biological path integration. This paper presents a scheme to classify path integration 
systems on the basis of the way the home vector records and updates the spatial 
relationship between the animal and its home location. Four extended classes of 
coordinate systems are used to unify and review both canonical and neural network 
models of path integration, from the arthropod and mammalian literature. This 
scheme demonstrates analytical equivalence between models which may otherwise 
appear unrelated, and distinguishes between models which may superficially appear 
similar. A thorough analysis is carried out of the equational forms of important 
facets of path integration including updating, steering, searching and systematic 
errors, using each of the four coordinate systems. The type of available directional 
cue, namely allothetic or idiothetic, is also considered. It is shown that on balance, 
the class of home vectors which includes the geocentric Cartesian coordinate system, 
appears to be the most robust for biological systems. A key conclusion is that 
deducing computational structure from behavioural data alone will be difficult or 
impossible, at least in the absence of an analysis of random errors. Consequently it is 
likely that further theoretical insights into path integration will require an in-depth 
study of the effect of noise on the four classes of home vectors. 

Keywords: navigation geocentric egocentric allothetic idiothetic 

Introduction 



Path integration (PI) flMittelstaedt and MittelstaedtUl980l : Mittelstaedt . 1983b) 



(IMoller and Gornerl . Il994l ) , birds ( Ivon Saint Paull . Il982l ) , crabs (ILayne et all l2003a! Jb 



is a navigation strategy many animals are capable of using, including ants (|Cheng et al. 
20091: iMuller and Wehnerl . Il988l iRonacherl. l2008h . bees (Ivon Frischl. Il967f). spider s 
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Zeil, 



1998), rodents fjMittelstaedt and Mittelstaedtl . ll980[ ) and humans fjMittelstaedt and Glasauerl . 



199ll ). whereby the animal maintains an estimate of its location as it moves around, 
by integrating its velocity over time. The animal's estimate of its location is re- 
ferred to as the home vector (HV), since it can be thought of as a vector connecting 
the animal's current location to the starting point of its journey, and allows it to 
return directly to the sta r ting poin t . For reviews i ntrod u cing; this behaviour see 



Collett and Collettl feOOOf ); iGallistell (11990f ): iRedishl (119991 ); IWehner and Srinivasan 
fl2003h . 

This paper undertakes a systematic comparison of several alternative ways of de- 
scribing PI, based on the coordinate system in which the HV is expressed to define 
the spatial relationship between the animal and its home location. Previous math- 
ematical models of PI have generally chosen a single coordinate system or frame 
of reference, often on the basis of incomplete or implicit assumptions. The only 
previous attempt to tran s late m odels between the alternative coordinate systems 
(IBenhamou and Seguinotl . Il995l ) made a weakly justified assumption that PI should 
only be thought of using one specific reference frame. To complicate matters, most 
published neural network models of PI were not defined using the standard coor- 
dinate nomenclature. This has resulted in a wide range of seemingly inconsistent, 
even contradictory opinions concerning the necessary properties of a biological PI 
system. At various times, some version of each of the standard coordinate systems 
have been proposed to be the foundation of the correct model of PI in the arthropod 
literature. 

In systematically comparing multiple coordinate systems, we explore and at- 
tempt to resolve several important questions. For example, the HV, being a vector, 
is suggestive of a polar representation. On the other hand, once near home, po- 
lar representations have significant drawbacks. An egocentric representation seems 
most intuitive for heading home. However, the egocentric position of home moves 
with every turn, even without displacement. Are these types of issues peculiar to 
specific coordinate systems? Do they persist following coordinate transforms? Can 
all major aspects of PI be represented equivalently in all coordinates systems? Are 
there theoretical clues as to which might be favoured biologically? It is important to 
recognise at the outset that there are intervening steps between the maintenance of 
a HV and its utilisation for finding home. For example, the HV could theoretically 
be represented in a geocentric framework (for definitions see below) but the motor 
commands needed to follow that HV may be more suitably transformed into an 
egocentric framework. This and related issues will be discussed later. 

This paper covers the four standard coordinate systems within which PI has 
previously been modelled, and provides equations for translating HV values between 
them. It also introduces a scheme for classifying PI models into four extended 
families of coordinate systems for those models which do not easily fit into any of 
the four standard systems, allowing the classification of virtually any conceivable PI 
system. Based on the analytical results, a critique is given of the usage of particular 
coordinate systems in existing models. Equational models of biologically important 
aspects of PI are introduced in all four of the standard coordinate systems. These 
include HV updating, steering, searching, and systematic errors. This forms the 
basis for an objective analytical comparison of the properties of the coordinate 
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systems. This paper considers PI in the absence of noise, and lays the foundation 
for later work which will incorporate the effects of random errors on PI to compare 
the noise tolerance of different classes of HV. Specifically, the definitions of the four 
extended families of coordinate systems were chosen to be suitable both for the 
present paper and for the study of noise tolerance. 

The general conclusion reached is that, mathematically speaking, we confirm 
that all the coordinate systems can adequately and usefully describe PI (in the sense 
of documenting or giving insights into navigation behaviour). However, geocentric 
Cartesian-like systems (see below for definition) appear the most robust solution 
for implementing a full PI system (in the sense of modelling the way in which an 
animal's nervous system needs to process and update information), particularly 
when an allothetic compass is available. 



Classification of Coordinate Systems for Path Integration 

Before considering which coordinate systems are most suitable for describing or 
implementing PI, it is essential to be able to classify the system which a given model 
uses. This section begins with four well known, standard coordinate schemes, and 
generalises them into four extended classes which can be used to classify virtually 
any conceivable model which can carry out accurate PI. 

Most existing equational models of PI can be directly assigned to one of four 
'standard' coordinate systems according to the type of HV used. The class depends 
on whether the animal's position is given in Cartesian or polar form and whether 
the journey's starting point (the geocentric case) or the animal's body (the ego- 
centric case) is used as the origin of the system. This leads to the four standard 
coordinate systems: geocentric Cartesian (GC), geocentric polar (GP), egocentric 
Cartesian (EC) and egocentric polar (EP). This paper considers only PI on a flat 
two dimensional plane, hence the simplest complete HV contains two values. Such 
a system can be extended to include PI on a non-flat surface by taking account of 
the local gradient, without the need for a full three dimensional HV. Desert ants 
appear to use such an approach, wh ere an essen t ially t wo d imensional PI system 



is made to cope with uneven ground (jGrah et all 120051 . 120071 ). This paper will use 



the following symbols to express the four standard types of HV: GC as (x,y), GP 
as (r, 9), EC as (x',y r ) and EP as (r',9'). The symbol <p will be used to indicate 
the animal's compass heading measured anti-clockwise from the direction of the x 
axis. Fig. [T] shows the meaning of the four types of HV diagrammatically. Table 1 
gives the equations needed for converting any standard HV into any other. Table 2 
summarises all the abbreviations used in this paper. 

Geocentric coordinates express the animal's position relative to the ground, with 
the origin corresponding to the starting point of the outbound journey and the 
direction of the axes corresponding to fixed directions with respect to the ground. 
The geocentric frame of reference is a special case of an 'allocentric' or 'exocentric' 
one: alio- or exocentric coordinates are those defined relative to something external 
to the animal's body, but in this paper the external reference used is always assumed 
to be the ground (for example the model animal is never on a table-top which can be 
rotated relative to the ground), making geo- synonymous with alio- and exocentric. 
Egocentric coordinates express the home position relative to the animal's current 
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position and orientation. The origin in this case is the centre of the animal's body, 
the x' axis corresponds to the forward direction along its body axis and the y' axis 
to a perpendicular axis pointing to the left from the animal's point of view. 

This paper presents equational PI models from the four standard classes of co- 
ordinate systems, because they encapsulate the essence of the majority of existing 
models, and they are mathematically convenient for analysis. However , to pr operly 
classify all the neural network models reviewed below (except Jander ( 19571 )). it is 
necessary to broaden each of the standard coordinate systems into an 'extended' clas s 
of coordinate systems based on the vectorial scheme employed by ICheungl (120071 ) . 
The four extended classes are: geocentric static- vectorial (GS, the class containing 
GC); geocentric dynamic- vectorial (GD, containing GP); egocentric static- vectorial 
(ES, containing EC); and egocentric dynamic- vectorial (ED, containing EP) repre- 
sentations (abbreviations summarised in Table 2). The distinction between geo- and 
egocentric remains the same for the extended classes but we introduce the concept of 
static and dynamic vectors. Static vectors are a superset of Cartesian vectors, where 
the HV representation now contains two or more values which are all distances mea- 
sured in strictly fixed (static) directions within a geocentric or egocentric space. For 
example instead of measuring distances along x' and y' axes defined relative to the 
animal's body and separated by 90 degrees (a standard EC vector), we can measure 
three distances along axes defined relative to the animal's body and separated by 
60 degrees. This is then an egocentric static vectorial representation and a model 
using such a HV would be classed as ES. If instead one or more of the distances 
are measured along an axis whose direction (within the defined ego- or geocentric 
space) is dynamically updated with changes in heading, then the vector is dynamic. 
This is clearly the case for all polar coordinates, since the distance defined by the 
scalar coordinate, such as r, is in a direction defined by the angular coordinate, such 
as 9. It is clear that some key properties of the standard coordinate systems must 
also hold for their respective extended classes. A comprehensive demonstration of 
this is, however, beyond the scope of the present paper (but see Appendices C and 
D for a brief examination of two general classes of GS HV update equations). 

One further definition will now be introduced, relating to the nature of the direc- 
tional cues available to the PI system. While not part of the classification of the HV, 
the nature of these cues has important consequences for the information processing 
required by the PI system. In the present account directional cues will be classified 



as allothetic or idiothetic following the scheme presented in ICheung et al.l ( 12008 



Table 1), depending on how the animal maintains an estimate of absolute heading. 
If the animal is required to perform neural integration to obtain the estimate of 
heading, then the directional cue is idiothetic. If no integration is required and the 
cue indicates absolute heading directly then the cue is allothetic. For example, ro- 
tational optic flow is generated from 'allothetic' sensory data (the movement of the 
eye relative to objects in the external environment), but provides only a measure 
of rotation rate, which must be integrated internally to estimate absolute heading, 
and is therefore cla ssified here as an id iothetic directional cue. This definition differs 
from that used by iMittelstaedt (l2000h . who defined idiothetic as "spatial informa- 
tion that can only be gained by means of the agent's active or passive movement", 
thereby excluding optic flow, which could arise from the flow of objects in air or wa- 
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ter past a stationary animal. However, since rotational optic flow accumulates errors 
in an analo gous way to (idiothet ic) vestibular or proprioceptive cues, we adopt the 
definition of Cheung et al.l ( 20081 ). since it is a more convenient basis for the analysis 



of errors in PI, and still agrees with Mittelstaedt's definition in spirit. 

The geometrically correct methods for updating a HV in all four of the stan- 
dard coordinate systems are given in Table 3 (HV Updating column), presented as 
continuous-time differential equations. These equations show how to update geo- 
centric HVs using allothetic directional cues and egocentric HVs using idiothetic 
directional cues, since these are the most concise forms. For comparison and com- 
pleteness, Tables 4 and 5 show the continuous-time equations required to update 
HVs using the alternative kind of compass input, and also present geometrically 
correct equations in discrete-time for all cases. 

Coordinate Systems in Existing Models 

The choice of coordinate system must be of fundamental importance for un- 
derstanding animal PI, since it determines the neural architecture, the nature of 
the information processing, and perhaps the type of sensory input required to per- 
form Pi-related behaviours. However, it is rare to find a complete discussion of the 
subject in published PI models. Typically, a published model adopts one partic- 
ular coordinate system without reviewing the alternatives. Therefore, one might 
consider many of the existing models to be proofs of concept or existence proofs, 
where such a review of alternatives could be considered as beyond the scope of the 
study. One obvious difficulty in the arthropod literature is the lack of neurophys- 
iological data, particularly in behaving animals. The paucity of data does allow 
room for speculative modelling. However, an existence proof cannot tell us whether 
behavioural consequences are specific to the assumed PI coordinate system. In 
contrast, the simple underlying geometrical structure of PI (limited admittedly by 
its noisy and approximate biological implementation) leads to simple methods for 
translating from one coordinate system to another (see Comparison of Coordinate 
Systems). Hence we should expect to be able to think about and describe a PI 
system from different reference frames, even if its actual neural implementation is 
in one specific system. Furthermore, a consideration of PI in alternative reference 
fra mes may lead to insights not obvious in the origin al coordinate system. 

Jander ( 1957 . in German, summarised in English in Maurer and Seguinot ( 1995h ) 



introduced the first published mathematical model of PI. The model assumed that 
the animal moved at a constant speed, making it simpler than most subsequent 
models, which generally provided some means to accommodate variable speeds. It 
also contained no true HV. The animal instead maintained a time-weighted average 
of its heading relative to an external light source, giving a single scalar in place of 
the usual two valued HV. The scalar was effectively the angular coordinate of a GP 
HV. Such a choice was natural since, during any initial straight part of a foraging 
trip where the home location remains directly behind the animal, the angle to the 
light is equal to the angular part of a GP HV, but is unrelated to the angular part of 
an EP HV which is always 7r during this period. Unfortunate ly, although presented 



as being correct, Jander's model is geometrically erroneous (IMaurer and Seguinot 
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19951 ) and, although clearly related to GP PI systems, is not strictly a member of 



this class. 



T he Mittelstaedts' PI m odel, the bicomponent model, ( IMittelstaedt and Mittelstaedt 



19731 ; iMittelstaedtl . Il983bl ) used a GC HV. This choice seemed to stem from earlier 



work related t o orien tation, including the work of von Hoist on orientation in fish 
(IMittelstaedtl . ll983aT ) . Fish motor responses to deviations from an external orienta- 
tion cue were found to be approximately sinusoidal, thus providing feedback capable 
of stabilising the animal. A problem arose when considering how an animal may 
orient at some arbitrary angle to an exter nal cue. Mittelstae dt proposed his recipro- 
cal modulation of bicomponents theory (IMittelstaedtl . Il962l ) as the solution to this 
problem. The present and desired angular deviations from the external cue were 
each broken down into a sine and cosine component (hence 'bicomponent model'). 
The two pairs were 'reciprocally modulated' (the sine of one pair being multiplied 
with the cosine of the other pair and vice versa), then the values were summed, 
and used to control the rotation of the animal. For PI, the sine and cosine values 
from an allothetic compass cue (or sometimes an idiothetic compass generated from 
integrating rotation rate) were used as inputs. A GC HV was the natural choice, 
obtained by multiplying the compass inputs by the forward speed and integrating 
over time (see Table 3, GC HV Updating). Homing was performed using recipro- 
cal modulation, the system deriving its desired compass heading from the current 
HV (see Table 3, GC Homing). GC coordinates therefore stemmed indirectly from 
the observation of sinusoidal motor responses in animals, generalised using the re- 
ciprocal modulation hypothesis. Neurons showing actual sine and cosine response 
functions to the orientation of external stimuli have since been found in nature: in 



response to touch in leeches (ILewis and Kristan Jr.l . Il998l ) and air flow direction in 



crickets (IMiller et al.l . Il99l[ ) . The latter could in principal constitute an allothetic 
compass cue given a sufficiently stable wind direction. Skylight cues are generally 
considered to be the most important allothetic compass cue for insects, and here 
no simple sine wave responses are known. Instead, polarisation-sensitive (POL) 
neurons in the early stages of the polari sation vision pathway of c rickets give a si- 
nusoidal response with a period of 180° (IHenze and Labhartl . 120071 ) rather than the 
360° predicted by the bicomponent model. There is evidence that information from 
polarisation and spectral cues as well as the solar elevation are combined to pro- 
duce a reliable, time comp ensated compass signal, but as yet no explicit sine-wave 
signal has been detected (jPfeiffer and Hombergl . 120071 ) . A positive feature of the 
bicomponent model was its explicit inclusion, often lacking in later models, of all 
the information processing required to generate a HV from clearly defined compass 
inputs and to generate motor control signals from the HV capable of taking the 
animal back home. 



Miiller and Wehnerl (119881 ) used a PI model based on a GP HV, to explain PI in 



desert ants. This choice probably stemmed from the model's starting point (draw- 
ing inspiration from Jander's model) as a calculation of the arithmetic mean of the 
animal's compass heading weighted by the distance walked in each direction. The 
use of Jander's model as the starting point, rather than the Mittelstaedts' model, 
was motivated by the idea that the ant was performin g a computation that is math - 
ematically simpler than the geometrically correct one (IMaurer and Seguinotl . Il995l ). 
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In spite of its starting point, the final form of the model is simple to derive from the 
discrete-time approximation (see Table 4, GP Discrete) of the geometrically correct 
continuous-time GP equations (see Table 3, GP HV Updating) by using a piece-wise 
linear approximation of the cosine function, a polynomial approximation of the sine 
function, and an additional scaling factor. These complications weakened the appar- 
ent link to Jander's GP-like model, and made the choice of coordinate system seem 
more arbitrary. Observed desert ant PI behaviour suggests a potentially disastrous 
drawback of using polar coordinates. The ant's HV is assumed to be set to zero 
magnitude (r = 0) when the animal initially leaves home to begin foraging. During 
homing the animal is assumed to be aiming to reach the goal location indicated by 
the HV (the location where its HV will have reached zero m agnitude again) , and dur- 
ing se arch behaviour it repeatedly passes over this location fjWehner and Srinivasanl . 



19811 ). The polar equations produce an infinite rate of change in the angular HV 



component at this location, due to the division by r (or r' for the egocentric equa- 
tions, see Table 3, GP and EP HV Updating) and more generally, the closer to 
home the animal gets the larger the c hanges in the angular compo nent must be- 
come. The Hartmann-Wehner model ( iHartmann and Wehnerl . Il995[ ) was a neural 
network implementation based on polar coordinates (see below). One important ad- 
vantage of the Muller-Wehner model was that its parameters could be adjusted to 
reproduce the observed systematic deviations from direct homing when desert ants 
were forced to walk along L-shaped channels ( Miiller and Wehner . 1988) (see below 
for a more detailed discussion of this phenomenon). However, it remains uncertain 
how generally this behaviour relates to natural navigation. Furthermore, as alluded 
to previously, it is not clear whether this ability of the model supports the argument 
for a GP HV for biological PI since it was not shown whether a translation of the 
model into other coordinate systems would alter its ability to reprod uce the errors . 



Finally, this model is not the on l y one able to account for the errors (IMerkle et al. 



2006M IVickerstaff and Pi Paolo! . 12005b . 



Gallistell ( 119901 ) presented two PI models in the form of equations and flowchart 



diagrams. As a prelude, Gallistel discussed the need to choose between polar and 
Cartesian coordinates for modelling PI, both apparently of the geocentric kind 
("Does the animal represent its position relative to the nest on a rectangular grid 
or as a certain distance and angle away from the nest?" p. 71). Gallistel considered 
polar coordinates to be appealing because they allow the HV to be stored in a form 
closer to that required for homing, but that they were less computationally attractive 
due to the presence of positive feedback loops in the updating system (although no 
quantitative treatment of this claimed weakness was given). The discussion lacked 
any explicit reference to the choice between geocentric and egocentric coordinates. 
The first of the two PI models presented used GC coordinates and was similar to the 
Mittelstaedt model, except that it used a more complicated homing system. The 
assumption appeared to be that the PI system must output the distance and turning 
angle required to reach home from the current location. To find these the system 
used a square root function and additional angular calculations. The Mittelstaedt 
model had already clearly shown how to accomplish homing without these extra 
complications, using the reciprocal modulation method (see Table 3, GC Homing). 
The second of the models presented was a polar PI system. From the definition 
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of the symbols used, the HV was intended to be of the EP kind, but the update 
equat ions given were not correct for this kind of HV (as iBenhamou and Seguinot . 



19951 . have observed). To be correct GP update equations, both right-hand sides 



wo uld need to be multiplie d by negative one. 



Benhamou et al.l ([1990) discussed the effects of noise on PI, which was referred to 



as 'the egocentric coding process'. The authors concluded that PI systems are much 
more sensitive to noise from idiothetic directional cues than either noise from allo- 
thetic compass cues or noise from estimating forward speed. The introduction set 
out a dichotomy between exocentric and egocentric coding processes. An exocentric 
coding process meant landmark-based navigation, where goal locations were mem- 
orised relative to landmarks. An egocentric coding process consisted of processing 
'route-based' (velocity cue) information, and was stated to be largely synonymous 
with the terms dead-reckoning, inertial navigation, path integration and route-based 
navigation as used by other authors. The authors were therefore distinguishing lo- 
cations memorised in geocentric and egocentric terms. PI was assumed to function 
using an EP representation of one or more goals (hence 'egocentric coding process'), 
justified by the supposed lack of dependence between individual stored goal coordi- 
nates in this system. The idea was to avoid positing a map-like representation of the 
animal's entire home range. However, if two or more EP goal coordinates were up- 
dated using the same velocity cues, they must lie within the same coordinate system, 
and therefore be dependent in this sense. Yet, given that egocentric coordinates must 
be constantly updated as an animal moves, the relative positions of two landmarks 
may diverge through cumulative errors thereby reducing their spatial relationship 
in representational space. This is certainly not a positive feature of EP coordinates 
for PI. Conversel y, geocentric (exocent r ic) co ordinates do not automatically lead to 
a cognitive map. ICollett and Graham! ( 12004] ) argue that if an animal can associate 
PI coordinates to an array of visually defined places, then by this fact it possesses 
a cognitive map, which Benhamou et al. wish to avoi d. Evidence f r om d esert ants 
suggests they cannot perform this kind of association ( ICollett et al.l . 120031 ). but this 
clearly doesn't invalidate the many geocentric PI models for this animal, which do 
not rely on the association of HVs with landmarks. In conclusion, the justification 
given by Benhamou et al. for using egocentric coordinates was weak, whereas the 
ch oice of polar coordinates was mad e without any clear justification. 



Wittmann and Schweglerl (119951 ) presented a neural network model of PI where 



the HV was stored on an array of neurons called a sinusoidal array. This design 
originated from hypotheses about how rats might encode the distance and direction 
to landmarks, but had the distinct advantage of making vectorial addition rather 
trivial, and it was also noise tolerant. The activity pattern along the array described 
a sine wave whose wavelength spanned exactly 360 degrees of angular space, and 
whose amplitude and phase encoded the distance and angle of a geocentric HV 
respectively. At first glance, this model resembles a polar representation. However, 
for present purposes the model must be classified as GS (i.e. a superset of GC). 
To understand this, we need to consider the HV updating process. Each step was 
represented as a sinusoidal activity pattern, and the net HV was generated from the 
cumulative sum of the activity patterns of all steps taken. Each neuron in the array 
represents a fixed geocentric compass direction. Any two neurons separated by one 
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quarter of the length of the sinusoidal array could be interpreted as GC coordinates 
on a certain set of fixed geocentric axes. The activity level and location of the 
most excited neurons in the array indicated the linear and angular component of 
the HV in GP form, witho ut having to per form anything resembling a square root 
or arctangent function, (cf. iGallistell . I1990I ). Hence the information in the HV was 
essentially available simultaneously in Cartesian and polar forms. The use of a 
geocentric rather than an egocentric HV arose naturally from the use of an array of 
input neurons representing an allothetic compass cue. The neurons of the sinusoidal 
array HV then essentially accumulated input from the corresponding geocentric 
compass input directions. The apparent contradiction of a GP HV readout from a 
GS representation stems from using redundant properties within the one PI system. 
In a GS representation, accurate PI can occur in at least two different ways. For a 
small number of static vectors, accurate PI requires a particular transformation of 
the heading direction prior to updating the HV (see Appendix C). However, for a 
large number of static ve c tors, the transformation is irrelevant (see Appendix D). 
Wittmann and Schweg ( 119951 ) used a large array of neurons and also assumed a 



sinusoidal activity pattern simultaneously, thereby carrying redundant information 
about the HV. 



Hartmann and Wehnerl (119951 ) presented a neural network model of PI which set 



out to reproduce the systematic homing errors seen i n desert ants in specific exp er- 
imental conditions, taking the Miiller-Wehner model (IMuller and Wehnerl . Il988l ) as 
a starting point. Howe yer, instead of deriving the model from a distance- weighted 
mean-heading model as IMuller and Wehnerl (119881) did, it was introduced as a mod- 
ification of geometrically correct GP PI (see lHartmann and Wehnerl . Il995l . Eqns 4 



and 5). Although this foundation contained minor inaccuracies (first, the paper dealt 
exclusively with angular units of degrees but Eqn. 5 was correct only for radians; 
second, the equations given were in fact only approximations to the geometrically 
correct discrete-time PI equations, see Table 4), these would probably not have pre- 
vented the model from performing PI quite accurately under most conditions. The 
neural implementation of these equations was then made deliberately imprecise, in a 
manner which mostly mimicked the mathematical approximations employed by the 
Mi iller-Wehner model, i n orde r to reproduce the systematic homing errors reported 
by IMuller and Wehnerl (119881 ) . In the neural network of lHartmann and Wehner 
(119951 ). the HV was explicitly separated into magnitude and angular components in 
the manner of polar coordinates. The distance magnitude was in the direction indi- 
cated by the variable angular component, hence the model must be classified as GD, 
a clear consequence of being based on the GP Miiller-Wehner equational model. The 
HV magnitude was represented by how far activation reached along a linear chain of 
neurons, and the HV direction was represented by the position of a patch of activity 
on a circular chain of neurons. Unlike the Wittmann-Schwegler model, informa- 
tion related to the Cartesian translation of the HV was not simultaneously available 
within the network using this representation. This direct implementation of a polar- 
like HV appears to inherit the divide-by-zero feature of the Miiller-Wehner model, 
requiring an infinitely fast rate of change of the angular HV component at the nest. 
The authors assumed that this value would change only very slow ly "as long as 
the ant does not run around its nest" (lHartmann and Wehnerl . Il995l . p. 488), but in 
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fact the problem would be if the ant ran over its estimat e of the nest location, a 
beha viour which is observed during ant search behaviour (IWehner and Srinivasanl . 
198ll ). The activity patch representing the HV angular component would be re- 



quired to jump discontinuously by a hal f turn, whi c h see ms impossible given the 
shifter circuitry employed by the model. I Chapman! (119981 ) reported problems get- 
ting the model to function correctly during tests on a mobile robot, particularly 
during homing, but did not state that this problem was limited to the vicinity of 
the nest, suggesting perhaps a further but as yet poorly defined problem of this 
model. One possibility may be the fact that approximate update equations can lead 
to negative values of magnitude which cannot be sustained by the aforementioned 
neural chain architecture. 



Benhamou and Seguinotl (119951 ) studied five PI models with the aim of trans- 



lating them into discrete time recurrent EP formulas. It is clear that the authors 
considered the EP reference frame to be the most appropriate for models of PI, but 
this choice was not adequately justified. The rationale given was mostly tautological: 
"From the animal's point of view, memorising the home location... is likely to be 
an egocentric coding process" (p. 463). Two references were given to papers which 
suggest that the caudate nucleus may be involved in processing egocentric positional 
information. This is weak justification since this brain area may not be involved in 
PI (or at least may not be the core of PI), and in any case three of the five models 
considered were originally formulated for arthropods which lack this brain area all 
together. Despite the weak basis for choosing the EP coordinate system, the authors 



felt justified in dismissing the bicomponent model ( IMittelstaedt and Mittelstaedt 



1973 iMittelstaedtl . Il983bf ) on the grounds that it cannot be easily converted into 



EP according t o thei r methodology. 



Benhamoul (119971 ) considered systems of reference for navigation in general, in- 



cluding PI. He noted that geo- and egocentric formulations of PI are in some sense 
mathematically equivalent, but claimed that egocentric formulations do not require 
an external reference direction whereas geocentric ones do. This is incorrect. It is 
possible to store the HV in geocentric terms but t o employ idiothetic dire c tiona l 
cues. For example the Hartmann-Wehner model ( Hartmann and Wehnerl . 1995) 



used this approach, as did a version of the bicomponent model (IMittelstaedt 



19851). 

An egocentric HV can also use information from an allothetic compass. Tables 4 
and 5 give details of how either ki nd of directi o nal cu e can be used with either kind 
of PI system. Even if the claim of iBenhamoul ( 119971 ) was true, it could not provide 
much justification for using an egocentric reference frame since ICheung et al.l ( 120071 . 



20081 ) have shown clearly that an external compass leads to more accurate naviga- 



tion than using purely idiothetic cues. Benhamou stated that when considering the 
neurobiological basis of PI "it seems better to adopt the animal's point of view, and 
hence to consider that path integration in v olves an egocentric localisation system" 



(IBenhamoul . Il997l . p. 151). IWehner et al.l (Il996h is cited in s upport of this view. 



despit e the fact that the latter refers to the model suggested by iMiiller and Wehner 
(|1988|). who prese nted a geocentric PI model. 



Maurerl (119981 ) presented a neural network PI model employing an EP-like HV. 



The magnitude and directional components of the HV were stored in separate neural 
structures, hence the HV can be classified as ED (since it contained more than 
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two values representing the polar HV, it was not strictly a simple EP HV). The 
magnitude was indicated by the position of an activity packet on a linear array of 
neurons and the direction by the position of an activity packet on a circular array 
of neurons. The network was trained using back propagation and had no homing 
m echanism - the output of the model was simply the HV itself. Unlike the model 
of lHartmann and Wehnerl (119951 ) . the angular component of the HV was not moved 
by dedicated shifter circuitry, and may have been able to meet the need for this 
part of the HV to jump instantaneously by a half turn were the animal to pass over 
the estimated home location (although this ability was not tested). This network 
model was compared with one which (rather like Jander's model) did not contain 
a representation of the HV magnitude, consisting of only the angular component of 
the first model. 



Kim and Hallaml (120001 ) outlined a neural-network based PI system. The HV was 



stored on a circular array of neurons, each corresponding to a fixed compass heading 
(or rather to a small contiguous range of compass headings), hence the model was 
GS. When the animal was travelling in a direction corresponding to a given array 
neuron, it was assumed to accumulate activity in proportion to the distance the 
animal travelled in that direction. The mechanism mapping directional and speed 
inputs into this integrative behaviour of the array neurons was not specified. Unlike 
the Wittmann-Schwegler model a neuron's activation level did not decrease when 
the animal was heading opposite to its preferred direction, meaning that the state of 
the HV neurons actually retained some information about the path taken to reach a 
given location, in addition to the location itself. Also, since the activity pattern was 
not sinusoidal, there was no redundant GP information in the HV representation 
(see overview of Wittmann-Schwegler model above). Making the sam e apparently 
unnecessary assumption (see Table 3 GC Homing) as iGallistell (119901 ) the authors 
indicated that they assume the HV must have a read-out mechanism providing the 
distance and direction to home in polar form. In line with this a second circular 
array of neurons was proposed which processed the HV stored in the first array to 
yield a single active neuron pointing in the direction of home, with the magnitude 
of activation indicating the distance. The way in which this signal was to control 
motor output was specified only using a fragment of pseudocode, rather than neural 
mechanisms. 



Bieglerl (l2000[ ) classified PI systems with respect to their level of capability, such 



as whether they could store more than one goal location, and how landmark-based 
positional fixes could be used to increase the accuracy of navigation. Coordinate 
systems were covered briefly, including a classification of HVs into polar, Cartesian 
and map-like, and as egocentric or exocentric. Biegler's classification of map-like 
HVs as a class by itself contrasts with the approach taken here, where they will be 
classified as GS. This choice is motivated by the fact that fixed directions within 
the abstract map-like sheet of cells in map-based PI models correspond to fixed 
directions in the geocentric reference frame. In general agreement with the present 
account, Biegler considered all types of HV representation as strictly mathematically 
eq uivalent. 



Vickerstaff and Di Paolol (120051 ) produced a neural network model of PI usin g 



a genetic algorithm (a system for producing artificial evolution (iGoldbergi . Il989l )). 
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The best evolved network obtained was found to be closely related to the Mittel- 
staedt PI model when analysed, and used a GC HV. This was the only type of HV 
encountered in the various runs of the genetic algorithm. The failure of evolution 
to utilise other forms of HV probably stemmed from the choice of the compass in- 
puts made available to the networks, which were either of the exact kind or similar 
to those required for the GC HV updating equations of the Mittelstaedts' model 
in most of the experiments. Unfortunately, there were no experiments which used 
rotation rate sensors in place of compass sensors. The equational form of the ego- 
centric HV update process appear simplest when rotation rate <f> is the directional 
cue, while geocentric update equations appear simplest when an absolute compass 
heading <ft is available (see Table 3, HV Updating). Hence one may not have ex- 
pected an egocentric HV to be easily evolved under the evolutionary circumstances 
modelled. In detail, there may also be some doubt concerning the biological plau- 
sibility of the evolved neural network due to the use of a non-standard neuron 
feature which allowed the network to perform multip licati ons more easily. However, 
for co mparison, both lHartmann and Wehnerl (119951 ) and IWittmann and Schwegler 
(119951 ) included an ad hoc multiplication-like mechanism in their network models to 
accommodate the effects of variations in the animal's speed into the HV updating 
process. The neuron- model (IBeer and Gallagher! . Il992l ) used in the neural network 
model of IVickerstaff and Di Paolol (120051 ) was based on a leaky integrator equation, 
so that the HV representation decayed continuously. Due to the artificially imposed 
selection of efficient homing, the network evolved to compensate for the decay of the 
HV so that no systematic navigation errors occurred. When this compensation pro- 
cess was disabled manually, the resulting systematic navigation errors were found to 
closely resemble those report ed in desert ants und e r the specific experimental condi- 
tions tested and reported by iMiiller and Wehnerl (119881 ). The network also evolved 
to perform a kind of searching behaviour when the simulated agent had returned 
to the vicinity of, but had not yet found, the nest. The features of a HV stored 
on leaky integrators, and of the searching behaviour displayed by the model will be 
considered in detail in the next section. 



Merkle et ajj (l2006bf ) presented an equational PI model using an EC HV. The 



authors first reviewed existing models with respect to their choice of coordinate 
system, and noted that EC had previously been neglected. The simple form of the 
EC HV updating and homing equations (Table 3) were seen as an advantage, as 
was the use of egocentric coordinates, because the animal "perceives all sensory 
inputs relative to its own position and orientation" (p. 388, ibid). While this is 
certainly true, it does not necessarily follow that egocentric coordinates are the 
natural choice for PI. According to the geometrically derived PI equations considered 
in the present paper (Table 3, HV Updating) an allothetic compass cue (denoted <f> 
in the equations) giving the animal its absolute heading, is more directly relevant 
for updating a geocentric than an egocentric HV. The rate of change of the absolute 
heading would have to be calculated in order to update an egocentric HV (see Table 
5), implying a more complex computation. The equations show that the perception 
of its own forward speed (s in the equations) is equally relevant to ego- and geocentric 
PI. The perception of rotation rate (0) appears more relevant to ego- than geocentric 
PI, since the rate must be integrated to estimate absolute heading before being 
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used to update a geocentric HV (see Table 4), but has been shown to be a much 



less reliable cue than an allothetic compass (iBenhamou et al.l . Il990l ; ICheung et al 



20071 . 120081 ). A further potential disadvantage of using EC coordinates is that the 



HV must be updated in response to rotation of the animal as well as translation. 
The further the animal is from the home location the greater the magnitude of 
the changes caused by rotation: a full turn of the animal causes the HV value to 
perform a complete circle in representational space. If rapid changes to the HV have 
an energetic cost or are associated with an increase in errors then an egocentric HV 
would be disadvantaged for this reason alone. This line of argument predicts that 
an animal using an egocentric HV might find it advantageous to limit it rotations, 
especially when far from the nest, but that one using a geocentric HV would not 
have any such incentive on this account. We save a more detailed treatment of this 
line of argument for future work, but note that desert ants certainly follow t o rtuou s 
paths during foraging;, i nclud ing when they are far from the nest (jWehnerl . 120031 ). 



Wehner and Srinivasanl ( 119811 ) also found a high maximum rotation rate (during nest 



search) in Cataglyphis albicans ants of 4000° per second. Merkle et al. undertook 
an investigation of a number of possible explanations of systematic homing errors 
in desert ants, and concluded that (amongst other candidate models) the use of 
leaky integration with their EC HV could explain the data well. As the present 
paper will show, the use of leaky integration in any of the four standard coordinate 
systems produces exactly the same error patterns. Therefore this fact does not lend 
much support to the hypothesis that EC coordinates are being used by these species. 
Other plausible explanations of PI systematic errors may yet arise and it remains to 
be seen whether they may be equivalently represented in all four coordinate systems. 



Haferlach et al.l (120071 ) presented a genetic-algorithm evolved neural network 



model of PI, which used a GS HV. They provided the netw ork with the same type of 
compa ss sensors as the Vickerstaff-DiPaolo model. Unlike IVickerstaff and Di Paolo 
(120051 ). the network topology was constrained during evolution to force each compass 
sensor to have its own directly attached integrative neuron. Both of these features 
probably helped predispose the system to find a GS solution - the genetic algorithm 
found the simplest solution, namely each integrative neuron integrated the output 
of its compass sensor. Since their simulated agent travelled at a fixed speed, this 
is sufficient to provide a GS HV (i.e. no input from a speed sensor is required). As 
the authors note the network did not store magnitude and angular HV information 
separately, and is only superficially similar to the Hartmann- Wehner model in this 
sense. The main differences from the Vickerstaff-DiPaolo model with respect to 
HV updating are that the network was given three compass inputs instead of two, 
leading to a three part GS HV rather than a two part strictly GC HV, and like 
Jander's model, the simulated animal is assumed to move at a fixed speed. The 
homing mechanism also differed: whereas the Vickerstaff-DiP aolo model exploite d 



the 'reciprocal modulation' solution proposed by Mittelstaedt (IMittelstaedtl . Il962l ) 



the Haferlach et al. model used a method similar to that in the Wittmann-Schwegler 
and Hartmann- Wehner models, and seems to have avoided the need for an explicit 
multiplication operation in the homing mechanism. This suggests an hypothesis 
about the problem Mittelstaedt was solving with reciprocal modulation: if the HV 
is expanded into a redundant encoding containing more than two values (GS rather 
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than the subset GC), a neural network can produce the correct homing direction us- 
ing a simple repulsion-like mechanism (where the animal's heading is repelled from 
the direction of maximum HV magnitude) without needing to implement recipro- 
cal modulation (i.e. multiplication). Vickerstaff-DiPaolo used a two value HV and 
appeared to require reciprocal mo dulation to generate d irect homing paths from 
all positions. In their first model, lHaferlach et al.l (120071 ) used a three-valued HV 
and appeared to generate direct homing trajectories without multiplication (direct 
homing was explicitly selected for during evolution through the genetic algorithm's 
fit ness function). 



Bernardet et al.l (120081 ) presented a neural network model of PI using a GS HV 



where a population of neurons was divided into 36 groups each representing a fixed 
compass direction. The number of active neurons in a group recorded how far the 
animal had travelled in this direction. Like the Kim-Hallam model, the activity of 
HV neurons could only increase during motion i.e., motion in the opposite direc- 
tion did not decrease activation levels. The model appears to be unique in that 
it included a mechanism for updating the HV for variable speeds of travel which 
did not require neurons to perform multiplication. Rather than having compass 
neurons with graded firing rates, compass neurons were either fully active or silent, 
so direction was represented by the location of a single active neuron in a circular 
chain. The speed was encoded in the firing rate of a 'velocity' (sic) neuron. A set 
of gater neurons passed the speed signal on to the HV memory neurons only for 
the memory group corresponding to the direction of the active com pass neur o n. No 
true homing mechani s m wa s presented. Instead the authors, like iGallistell ( 1990 ) 
and iKim and Hallaml ( 120001 ) , assumed that the H V memory neurons must provide 
a 'read out' mechanism that outputs the direction and magnitude of the HV. 

For co mpleteness, we briefly mention m ammalian PI models. The hippocampal 



place cells ( lO'Keefe and Dostrovskvl . Il97ll ) , medial ento rhinal grid cells (IHafting et al. 



20051 ) and the more diffuse head direction cell network ( jTaube et al.l . |l990a 



bf) stud- 



ied in rodents have been hypothesized to form t he basis of a PI system (Bura k and Fietd. 



2009; iBureess et all . 120071 : iMcNauehton et all . Il996l . 120061 : ISamsonovich and McNauehton 



997). Both the toroidal neural network model, reviewed recently by lMcNaughton et al 



(200 6[), a nd th e algorithmic oscillatory interference model, recently reviewed by 
Burgess! (120081 ) use map-like PI schemes in a geocentric reference frame. There 



is rapidly growing consensus that rodent PI HV updating occurs in the grid cell 
networks, which typically display geocentric firing patterns as rats moved about 
experimental environments. The rodent PI models differ from the arthropod PI 
models presented above in using an abstract two dimensional space, or chart, within 
which neural activity patterns formed one or more activity packets, the location and 
movements of which fulfilled the role of the HV. This obvious difference does not 
however put them beyond the scope of our classification scheme: using the extended 
coordinate system classification scheme, neural units which fire in response to being 
at specific positions in space are representing GS vectors. Map-based PI models 
are therefore a special case of GS PI where the static vectors have fixed length as 
well as geocentric direction, hence limiting their spatial receptive field to a point 
(or small region). A 2-D array of such neural units could theoretically cover the 
entire PI range with or without tessellation. Although not so far proposed in the 
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literature, we note here that an ES HV could also be implemented using a map-like 
scheme, where the activity packet position represented the nest's location relative 
to the current position and orientation of the animal. Thus map-based PI systems 
fit easily within the classification scheme introduced in the present account. 



Comparison of Coordinate-Systems 



This section considers several biologically i mportant aspects of PI (see Table 3): 



HV u pdating, systematic Pi-navigation errors ( iMuller and Wehner.ll988tlSommer and Wehneri . 
2004), Pi-mediated homing an d Pi-m ediated systematic searching (IMuller and Wehneri . 



19941 : I Wehner and Srinivasanl . Il98lh . These aspects of PI are used as the basis 
for comparing the different coordinate systems. Systematic errors and systematic 
searching are considered in some detail. Differential equations defining a continuous- 
time model of each of these aspects of PI (Table 3) are used as the main basis for 
discussion, for reference discrete-time versions of HV updating, homing and search- 
ing are also given (Tables 4, 5, 6). Each of the models is presented in all four of 
the standard coordinate systems, and in each case the four models were obtained 
through exact transformation of the GC version of the model into the other three 
reference frames, ensuring that the four are mathematically strictly equivalent. The 
four sets of equations for each coordinate system were also validated using numer- 
ical inte gration (Runge-K utta fourth order method, with a fixed time step of 0.01 
seconds; iDormandl . Il996l ) to ensure they generated equivalent trajectories under a 
range of initial conditions (data not shown). 



Updating the HV 

HV updating, which could also be referred to as "path integration proper", is 
the process whereby sensory cues relating to velocity are used to update the current 
HV (see Table 3, HV Updating column). The nature of the updating process will 
obviously depend on the nature of the HV and the available sensory cues. Nonethe- 
less, all types of HV updating require information about the animal's speed s and 
also some kind of directional information. The classification of directional cue as 
allothetic or idiothetic was discussed earlier. In the equational models considered 
here <fi refers to the animal's compass heading and is most naturally interpreted as 
an allothetic compass cue, whereas refers to the animal's rotation rate and the 
natural interpretation is as an idiothetic cue. The HV update equations listed in 
Table 3 were derived from the geometrically correct GC equations for HV updat- 
i ng (Table 3, GC HV Updating). These w ere taken from the bicomponent model 
(IMittelstaedt and Mittelstaedtl . Il973l 119821 ) . but since they are geometrically exact, 
they can be considered as mathematical definitions rather than modelling choices. 
These were then mapped into GP, EC and EP coordinate systems (see Appendix 
A), such that the four resulting sets of HV update equations were strictly mathe- 
matically equivalent, although it should be stressed that they differ with respect to 
their application (see below). 

From the HV update equations, it is simplest to express the updating of geo- 
centric HVs using information relating to the value of <fi and simplest to update 
egocentric HVs using the value of <ft (Table 3, HV Updating). Interestingly, this 
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is more than a purely mathematica l featu re (unlike e.g. the polynomial approxi- 



mations used by iMiiller and Wehnerl (119881 )). and appears to be a general property 



of PI with neurobiological implications: any HV in an egocentric reference frame 
requires information about rotation (rate), while a geocentric HV requires absolute 
heading information. This implies allothetic compass cues are more suitable inputs 
for a geocentric PI system, whereas idiothetic directional cues are more suited to 
egocentric PI systems. However, this reasoning is not absolute since it is also possi- 
ble to use an allothetic compass to measure rotations and maintain an egocentric HV 
(Table 5, Allothetic directional cue column), or to integrate/sum idiothetic rotations 
to estimate absolute heading for maintaining a geocentric HV (Table 4, Idiothetic 
directional cue column). Nonetheless, the aforementioned directional-input/HV- 
coordinate combination seems the most mathematically parsimonious. 

Another feature which differentiates the mathematical HV update equations is 
the variation in the required rate of change in the HV as the animal moves around. 
As previously noted, the division by r or r' in the polar update equations (Table 
3, HV Updating) means that the rate of change of the angular coordinate increases 
towards infinity the closer the animal gets to home. This is a general feature of polar- 
like dynamic vectorial HVs. By definition, dynamic vectors have variable angular 
components which need to be updated during PI. At zero vectorial length, a dynamic 
vectorial angle is undefined. In contrast, the EC HV update equations show an 
increase in required rate of change in response to rotation the further the animal is 
away from home. This is a feature of ES HVs in general. By definition, the egocentric 
HV is expressed relative to the animal's body axes and therefore must dynamically 
account for every rotation of the animal itself (unlike geocentric systems) . Since the 
vectors are static, rotation implies remapping of current amplitudes. The longer the 
HV, the larger the ES vectorial amplitudes, thus requiring a greater absolute change 
in amplitudes necessary for a given rotation. In fact, only the GC equations have 
a maximum rate of change, proportional to the animal's maximum forward speed. 
This is true of GS PI systems generally (as can be see for the two general classes of 
GS HV update models shown in Appendices C and D). If rapid changes in the HV 
are associated with increased error accumulation or energetic costs, GS HVs should 
be favoured by biological PI systems. 

A distinguishing characteristic of GC/GS PI systems is that the HV requires 
only information from speed and directional cues, whereas the other three systems 
require information about the current state of the HV itself (e.g. for EC HVs, to 
update the x' coordinate requires information about the y' coordinate as well as 
speed and rotation rate). This suggests that in non-GS PI systems, HV information 
has to feedback via two pathways to give the baseline HV state, as well as the change 
of HV state. A GS PI system avoids this complication and should be considered an 
advantage. 

Reading the HV 

Some models of PI describe how to update the HV, but not how to use the HV to 
control the animal's behaviour. Others have a 'readout mechanism', where the HV 
is assumed to need translating into an alternative reference frame, but do not specify 
how this will then influence behaviour. This approach is in danger of deferring some 



16 



portion of the required information processing to an homunculus who looks at the 
readout and uses it to control behaviour. Here it is assumed instead that the HV 
must be used to directly control the animal's rotation rate and speed to perform the 
required Pi-related behaviour, the simplest of which is 'bee-line' homing. In this way 
all information processing required of the brain can be made explicit (although some 
Pi-related information process ing is beyond the scope of the current paper, such as 
ephemeris compensation, e.g. IWehner and Miillerl (1199310 . The homing equations 
used here were taken from the GC Mittelstaedt model (IMittelstaedtl . 120001 ) (see Table 
3, GC Homing), and were then mapped into all four standard coordinate systems 
(Table 3, Homing column) (see Appendix A). This homing model can be seen as 
more of an arbitrary modelling choice than the HV update equations, since there 
is no single geometrically correct way to control motion to produce homing. The 
Mittelstaedt model is arguably the simplest possible, and is easiest to understand 
in EC form (although originally given in GC): in EC the turning rate equation may 
be read as 'turn left while home is left, turn right while home is right' and the speed 
equation may be read 'go forwards while home is in front, go backwards while home is 
behind'. The equations will be considered here as generic steering equations, rather 
than a detailed biologically realistic model. The equations show that egocentric 
HVs contain information in a form directly usable for homing (for this simple model 
at least, and probably for homing-related behaviour in general), whereas geocentric 
HVs need additional information from the compass and are effectively dynamically 
translated into egocentric terms before being used for homing. 



Leaky Integrator Model of Desert Ant Systematic Navigation Errors 

Several studies have addressed the phenomenon of systematic navigation er- 
rors in animals navigating by PI, particularly in desert ants Cataqlyyhis forti s 
( Maurer and Seguinot . 1995 : Miiller and Wehner . 19881: ISommer and Wehner . 2004 ). 
In the experiments reported by iMiiller and Wehnerl (jl988l ) the ants travelled out- 
wards from the nest along a two-legged L-shaped channel to the feeder, and were 
then released to perform homing in a distant test field. In addition to the ex- 
pected scatter attributed to unbiased random errors, on average the ants homed 
in a direction different from the geometrically correct homing direction. The de- 
viations from the expected homing direction varied systematically with the angle 
of the outward L-shaped channels. The Miiller- Wehner model ibid, could explain 
these errors by assuming that the method used by the ant to update its HV dur- 
ing the outward portion of the journey was different from the geometrically correct 
meth od. Therefore this mode l , as well as subsequent attempts to explain these re- 
sults (jHartmann and Wehner . 1995: Merkle et all l2006bl ; IVickerstaff and Di Paolo . 
2005 : Iwittmann and SchweglerT 1995 ). attributed the homing behaviour to a kind 
of systematic error that the ants made when updating their HV. Interpretations of 
the Miiller- Wehner error should take account of the fact that they were observed in 
specific experimental settings and analysed in very specific ways, and thus may not 
relate to the ant's natural PI behaviour in a straightforward manner. In addition, 
to date, alternative explanations have not been rigorously excluded. Ant navigation 
is more than PI, and it is plausible that an optimal navigation strategy incorporates 
PI but does not follow its output exactly. For example the Australian desert ant 
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Melophorus bagoti trav els only half the h oming distance indicated by its PI system 
under some conditions ( INarendral . l2007af ). Optimal navigation may also not always 
entail travelling dir ectly towards the i ntended target, particula rly when additional 



cues are available (ISteck et all 120091 ; IWolf and Wehnerl . 120051 ). yet the PI system 



could still be geometrically exact per se. It is also possible that unbiased random 
errors at some level in an otherwise geometrically exact PI system may produce 
bia sed systematic error s in th e observed behaviour. 



Miiller and Wehnerl (119881 ) provided the first model of sy stematic errors in desert 



ants b ased on approximating the sine and cosine functions, lHartmann and Wehner 
(199 51) later produced a neura l network model based on the first model, while 



Wittmann and Schweglcr (199 5J) provided an alt ernative mechanism based on a sim- 



ple time delay, IVickerstaf f and Di Paolol ( 120051 ) applied a leaky integrator model 
using a GC HV and iMerkle et al.l (j2006bl ) undertook a systematic study of several 
candidate models including leaky integrators, using EC HVs. All of these models 
sought to explain the PI errors fro m the L-shaped chann els described above, with 



Vickerstaff and Di Paolol ( 120051 ) and IMerkle et al.l ( 12006bl ) indepe ndently proposing 



the sa me leaky integrator mechanism. In more recent experiments, ISommer and Wehner 
(120041 ) observed a different kind of systematic error in desert ants, that of underesti- 
mating the length of long straight jo urneys, and examined several can didate models 
again including a leaky integrator ( IMittelstaedt and Glasauerl. Il99ll. introdu ced a 



similar model for a human navigation task). IWittmann and Schweglerl ( 119951 ) pro- 
posed that underestimating the return distance, or at least beginning to search for 
the nest before running off the entire length of the HV, may have the selective 
advantage of increasing the likelihood of encountering familiar landmarks to guide 
homing. Recent work in Austra l ian des ert ants appears t o prov ide supporting evi- 
dence for these ideas ( INarendral . l2007al Jbl). IMerkle et al.l ( ]2006ai ) also show that, in 
Cataglyphis fortis, depending on the degree of the ant's uncertainty about the home 
location, systematic search behaviour does not necessarily start after the expected 
homing distance has been run off. 

This section uses leaky integration as a model (see also Appendix B) of systematic 
navigation errors in order to compare the four coordinate systems under conditions 
where HV updating is no longer geometrically correct. It will be introduced in 
GC form first (Table 3, GC Systematic Errors), but like the HV updating and 
homing models, mathematically equivalent equations for the three other standard 
coordinate systems are given (Table 3, remainder of Errors column). The GC form 
assumes that the geometrically correct sensory input values, s cos and s sin 0, are 
available as inputs, but that the integrator itself is 'leaky' and decays over time 
as a first order process with rate constant kn- Hence HV updating is no longer 
according to the geometrically correct formulae, and the animal could be considered 
as gradually forgetting earlier parts of its journey. Translation into EC (Table 3 EC 
Systematic Errors) reveals that, in order to exactly reproduce the behaviour of a 
GC PI system with leaky integrators, it is sufficient to make the EC HV coordinates 
behave as leaky integrators with exac tly the same decay rate. This result shows 
that the GC leaky integrator model of I Vickerstaff and Di Paolol ( 120051 ) and the EC 
model of IMerkle et al.l (l2006b[ ) are mathematically exactly equivalent (at least in the 
absence of noise), and predict precisely the same errors for all possible journeys. The 
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translations from GC and EC into GP and EP show that the decay can be considered 
to directly affect just the magnitude component of the HV since the decay constant 
only appears in the equations for r and r' and not for 8 or 6'. However, it should be 
emphasised that the HV of a non-straight journey can still end up pointing in the 
wrong direction, since a later step is effectively represented on a larger scale than 
an earlier step of the same physical length, leading to a distortion of the HV. 

It is perhaps surprising that deviating from the geometrically exact HV update 
equations in one coordinate system still yields concise, understandable equations in 
the other reference frames. In this case the leaky integrator concept maps easily 
into all four standard coordinate systems. This is probably not the case for all 
conceivable deviations from geometrical HV updating, but the present result serves 
as an existence proof that such equivalences do exist and may frustrate attempts to 
probe the inner work i ngs o f biological PI systems using only behavioural data (cf. 



Maurer and Seguinotl . Il995[ ) 



The leaky integrator model will now be fitted to two existing data sets from 
experiments with desert ants to determin e the decay rate paramet e r (thi s is a re- 
working of the same findings presented in IVickerstaff and Di Paolol . 120051 . using SI 
units). The model provides a good quantitative fit to the two sets of experimental 
data examined, although a different decay rate is req uired for each set. To mode l the 
systematic ho ming errors observed after L -shaped ( Miiller and Wehnerl . 1988 ) and 
long straight f jSommer and Wehnerl . 120041 ) foraging journeys it is sufficient to con- 
sider journeys made up of a small number of straight sections with uniform speed, 
where the differential equations have simple solutions without the need for numeri- 
cal integration. The simplest set of assumptions, which are nevertheless reasonably 
realistic for C. fortis ants (T. Merkle, personal communication), are that the an- 
imal walks with a constant speed (s) during foraging and homing, that it spends 
no significant time at the feeder before beginning to head home, that the decay 
rate is constant throughout the entire journey and that homing consists of walking 
straight towards the current estimated home location with no significant variations 
in heading or speed. The state of a GC HV after walking along a straight section 
with compass heading cf) at speed s for a time t will be: 



s cos <p \ , , s cos < 



k D J k D 

where the GC HV state at the start of the section is (xo,yo) and the decay rate is 
kjj. Upon reaching the end of the foraging part of the journey the ant will adopt 
a homing direction determined by the current HV state (x, y) such that it orients 
directly towards the estimated home location. Note this is a simple direct hom- 
ing path in a straight line, and is not the homing path predicted by the homing 
equations previously considered; the homing equations essentially show how the 
ant may adopt orientation ^homing without performing an explicit arctan calcula- 
tion (^homing = atan2(— y, —x) (where atan2 is the 4-quadrant arctangent function 
as defined in Table 1). Switching for convenience to the (strictly equivalent) EP 
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description, the value of the HV with the animal's orientation at ^homing must be 
r > = ^Jx 2 + y 2 , 9' = since home is now directly ahead. During homing only the 
value of r' will change according to the leaky integration model (Table 3 EP System- 
atic Errors), resulting in only one update equation which is simply r' = —s — k^r' 
because 9' = 0. The value of r' when homing first starts will overestimate the dis- 
tance the ant walks before the HV reaches zero due to the ongoing decay process. 
The homing distance will be: 



"fin 

k D 



k D r' + s 



where r' is the HV magnitude at the beginning of homing. This calculates the 
distance the animal will have walked before the HV magnitude reaches zero, which 
is taken as the perceived home location (see Appendix A for the derivations). 

A z ero decay rate yields geometri cally correct navigation, regardless of walking 
speed. lHartmann and Wehnerl (119951 ) showed figures of the observed homing direc- 
tions from three di f ferent L-shaped outward excursion experiments performed by 
Miiller and Wehner (1988) . An a verag e ant walking speed of 0.33ms -1 can be esti- 
mated from iMiiller and Wehnerl f 19881 . Fig. 1). Using this value for s, a decay rate 
constant of k^ = 0.0185s -1 is obtained by fitting the leaky integrator model to the 
homing direction data (data from the three experiments were pooled before fitting 
fco). Figs [2H show the predicted homing angles using this k D value. Assuming the 
same walking speed, a different decay rate constant {k D = 0. 00171s -1 ) is obtained 
by fitting the model to t he er rors observed during long straight journeys reported 



m iSommer and Wehnerl ( 120041 ) . Fig. [5] shows the predicted homing distances from 
this decay rate. 

Pendulum Model of Desert Ant Systematic Searching Behaviour 

When desert ants reach the end of their HV but do not find the nest entrance, 

This may be a 



they begin a searching behaviour (IWehner and Srinivasanl . 119811 ) 
method to compensate for the inevitable errors which accumulate in the HV dur- 
ing foraging excursions, and shows significant flexibility in response to the condi- 
tions of the preceding journey. Longer foraging journeys lead to decreased certainly 
about the nest positi on, which the ants compensate for by searching a wider area 
(IMerkle et al.l . l2006af). The spread of the search also increases the longe r an in divid- 
ual search lasts for (IMiiller and Wehnerl . Il994c IWehner and Srinivasanl . Il98l[ ). The 
ants can also learn to bias the shape of their search patter n in response repeated mis 



matches betw e en th e expected and actual nest position ( jCheng and Wehnerl . 12002 



Wehner et al.l . 120021 ) . Mathematically, the most efficient search pattern would be 
an ever expanding Archimedean spiral. However, if the ant failed to detect the en- 
trance on its first encounter, it would continue to spiral outwards and would never 
find it (IWehner and Srinivasanl . Il98ll ). Hence the real ants do not use this strat- 
egy, but perform a series of loops, gradually expanding overtime, but centred on 
and occasionally returning to the estimated home location. This means the amount 
of time spent searching per area decreases with distance from the centre, result- 
ing in a search density profile which is approximately symmetrical and bell-shaped. 
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Previous models of s earch i ng behaviour in arthropods inc lude iHoffmann (Il983al lbh: 



Miiller and Wehnerl (11994 ); IWehner and Srinivasanl (Il98ll ). 



To generate a model of searching behaviour suitable for a systematic comparison 
of the four coordinate systems, it would be convenient to have a simple modifica- 
tion of the homing equation s of th e Mittelstaedt model which enable it to produce 
searching. iKim and Hallaml ( 120001 ) have noted how their simple pseudocode-based 
homing system produced simple repetitive looping behaviour near the home loca- 
tion, and speculated tha t this might form the bas i s of a searching behaviour. The 
neural network model of IVickerstaff and Di Paold (120051 ) is very similar to the Mit- 
telstaedt model but also produced searching behaviour which is qualitatively similar 
to that of desert ants, including a radially symmetric bell-shaped search density. As 
discussed earlier, the model was produced using a genetic algorithm and was selected 
to perform Pi-mediated homing in the presence of noise, including sensory noise. 
This led to significant random errors in the HV. As a compensatory strategy the 
network evolved to perform a searching behaviour upon returning to the estimated 
home location. The full neural network model is too complex to use here but, for- 
tunately, analysis of the network (R. J.Vickerstaff, unpublished) produced a simple 
modification of the Mittelstaedt homing equations which generates searching in a 
manner virtually identical to that of the full neural network model. This section 
examines the model (Table 3, Searching column), in the original GC reference frame 
and also following translation into the other three standard PI coordinate systems. 

The model produces straight homing trajectories when the animal is far from 
home, but close to home a complicated non-repeating, looping search pattern de- 
velops (see Fig. E]). The search patterns still occur even though the noise used 
during the evolution of the model is re moved. Numer ical calculation of the maxi- 
mum Lyapunov exponent of the system (ISprottl . 120031 ) shows that the dynamics of 
the searching pattern are chaotic, indicating that the patterns do not have a simple 
repetitive structure, even in the absence of external noise. The search density in the 
presence of noise is a smooth radially symmetrical bell shape (data not shown), in 
the absence of noise the distribution is still roughly bell sha ped (see Fig. [71) . This is a 
simil ar search density profile to that observed in desert ants (IWehner and Srinivasanl . 
198 ll ). but the fine details of the trajectories do not appear to match those of the 



ants. 

While the original neural network operated using a GC HV, the searching dy- 
namics are perhaps easiest to understand in an egocentric reference frame. Since the 
search profile of the model is radially symmetrical around the home location (as it 
often is in the real ants) it makes sense to assume that the absolute compass heading 
at any point during the search is not important per se. The egocentric translations 
of the searching model (Table 3, EC and EP Searching) confirm this since they show 
analytically that the value of is not needed: the three variables (x',y',<p) form a 
self contained ('autonomous') system which encapsulates the complete dynamics of 
the search model expressed in egocentric terms, referring only to the rotation rate 
<f> and never to the absolute heading <p: 
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x' = <fiy' — s 

y = —<px 

= hy' - k 2 4> 



Fig. shows in EC-rotation space (i.e. (x',y',<p)) the search pattern shown in GC 
space in Fig. |6j Using parameter values k\ = 2.7973m~ 1 s~ 2 , k 2 = 1.308s -1 and s = 
lms -1 the model makes the animal spend half its search time within lm of the central 
nest location. Any such egocentric trajectory can be mapped into multiple geocentric 
searches since it is valid for any absolute orientation. To map an egocentric search 
trajectory to geocentric space requires choosing an initial value for the start of 
the time series and applying the appropriate conversion formulae from Table 1 to 
convert the HV to the corresponding geocentric values. Any value of can be used 
since the egocentric equations are independent of it. To produce the same search 
pattern using the geocentric equations requires numerical integration of not three but 
four variables: (x,y,<fi,(ft). In removing the dependence on absolute orientation the 
egocentric description employs one less variable but still encapsulates the complete 
dynamics and makes the underlying ordered structure of the search model more 
apparent. 

As previously noted the EC rotation rate control equation for homing was = 
k$y', described as 'when home is to the left turn left, when home is to the right 
turn right'. The equation for searching controls rotational acceleration rather than 
rotation rate directly and could be described verbally as 'when home is to the left, 
rotation tends towards the left, when home is to the right, rotation tends towards the 
right'. To gain a more detailed understanding of the searching model, consider the 
EP formulation of the same system (see Table 3, EP HV Updating and Searching): 



r' = — s cos 6' 



sin 9' 



k\r' sin 9' — k 2 c, 



Now consider the special case where forward speed is set to zero (s = 0), but the 
animal can still rotate: 



r > = o 

6> = -0 (1) 
4> = k ir 'sm9' - k 2 <p (2) 



Combining equations [T] and |2j 
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6' + k 2 6' + far' sin 0' = 



results in an equation exa ctly homologous to that describing a linearly damped 
pendulum f Strogatzl . 1994 ). with 'gravity' coming from the direction of home (as 
indicated by the HV) with a force proportional to the distance from home, r', and 
with k 2 as a linear damping constant. When not moving forwards the animal's 
orientation therefore acts as a damped pendulum which will always come to rest 
at a stable equilibrium pointing homewards, unless it is perfectly balanced in the 
direction pointing away from home, but this state is unstable. It is now easier 
to understand why the model produces direct homing when far from home, but 
searching behaviour close to home. Strong 'gravity' when far away quickly brings 
the orientation homeward, but when close the attraction is sufficiently weak to allow 
search loops to form. The searching pattern arises as the complicated interaction of 
the swings of the pendulum with the changes in the egocentrically perceived homing 
direction due to forward motio n. The values of fa and k 2 must be set correctly to 
generate search-like behaviour ( iVickerstafil . 120071 ) . 

The pendulum model's searching equations show that egocentric HVs contain 
all the necessary information to mediate searching, whereas geocentric HVs must be 
supplemented with compass information, whereby they are effectively dynamically 
translated into egocentric terms before being fed into the searching system. This is 
similar to the simpler homing equations. 

The pendulum model certainly seems to suggest that egocentric PI systems are 
well suited for mimicking the radially symmetrical Gaussian-like search profiles of 
desert ants observed under quasi-natural conditions. Whilst symmetrical searching 
may be simpler, perhaps even more intuitive to quantify in egocentric terms, it does 
not prove that the underlying PI system uses an egocentric HV. As shown in Ta- 
ble 3, symmetrical search-like behaviour can be generated regardless of the type of 
HV used. Furthermore, it has been demonstrated that ants may search for their 
nest asymmetrically if the ant repeatedly experience s an experimentally induced 



mismatch between its HV and the location of the nest f lWehner et all 120021 ) . Asym- 



metry develops in the searching distribution in the direction which compensates for 
the mismatch. The latter observation implies that searching behaviour is influenced 
by geocentric information. 

Nonetheless, the preceding analysis suggests that in practice, regardless of the in- 
ternal PI system, there are advantages in mapping searching and homing trajectories 
into an egocentric reference frame. This may assist in revealing underlying patterns 
in the data by removing the absolute orientation from consideration, particularly 
where the characteristic behaviour is independent of it, and to emphasise the steering 
control process. The simple pendulum model presented above also demonstrates the 
utility of chaotic dynamics in generating an efficient, non-repeating search pattern, 
and that such dynamics may arise spontaneously from simple deterministic rules 
using nothing more than the HV and current heading. Taken together, these results 
disprove the argument that complex searching must be a distinct behavioural pro- 
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gram to that of early homing. Furthermore, the presence of homing and/or searching 
behaviour per se cannot be used to support or oppose an argument for a particular 
coordinate system for PI. 

Storing Locations with the HV 

Navigation often requires a complex interplay of mechanisms of which PI is one. 
Although the mathematical analyses presented in this work do not relate directly to 
other forms of navigation, it is prudent to illustrate the potential relevance of the PI 
coordinate system to other forms of navigation. Here we briefly consider whether the 
use of landmarks in combination with PI for navigation may be more readily achieved 
in certain PI coordinate systems. There is evidence that rodents are able to associate 
landmark information with the state of the PI system, and use land mark informa- 



tion to correct errors by recalling the associated path integrator state (lEtienne et al. 



20041 1 . For this purpose, a geocentric map-like coordinate system possesses distinct 
computational advantages. Firstly, in a geocentric coordinate system, the coordi- 
nates of multiple landmarks may be stored in memory without the need for continual 
updating as the animal moved around. In contrast, if using egocentric coordinates, 
landmark coordinates must be updated even if the animal turned on the spot. Sec- 
ondly, in a map-like PI scheme, landmarks can become linked to the state of the 
PI system through associative learning in a path- independent m a nner. In contrast, 
i n som e GS HV representations, such as iBernardet et al.l (120081 ); iKim and Hallam 
(120001 ). the HV state is partially dependent on the path taken to reach the current 
location, because the activity levels of neurons are never actively reduced during 
locomotion. Thirdly, in a map-like scheme, a distinct subset of neural units may be 
associated to each landmark, reducing the possibil ity of ambiguity. In contra s t, eve n 
path-independent GS HV representations, such as lWittmann and Schweglerl (119951 ). 
require every neural unit to be at least partially active at every location in the PI 
range. Admittedly, polymodal firing fields of grid cells may partly compromise the 
advantage of the map-like GS scheme in this respect, since they introduce spatial 
ambiguity in the firing fields, each of which is active at multiple locations in space. 
Therefore binding a landmark to grid cell activity re quires connections to not just 



one but a population of grid cells ( IMoser et al.l . 120081 ) 



Although it is unclear in arthropods w hether the state of the PI system per 



se can associate with landmarks (but see ISrinivasan et al.l . 119971 ) . there is little 



(Collett et al.. 


1999; 


von Frisch 


. 19671 Rilev et al.. 


2005; 



Wolf and Wehneri . l2000h . 



If an animal stores a geocentric HV in memory while at a key location, the memory 
contains sufficient information for returning to the location by PI later. If an ego- 
centric HV is used for this purpose, the animal must also store its compass heading 
when at the location, otherwise the stored HV does not uniquely specify a geocen- 
tric location, but rather a circle of some radius centred on home. Alternatively, 
the animal must maintain a separate egocentric HV (or rather goal vector) for each 
goal location and constantly update these while moving. A more general property is 
seen in Table 1, where transformations between egocentric to/from geocentric refer- 
ence frames requires the heading 0. In contrast, transforms between static to/from 
dynamic vectorial representations do not require 0. Therefore, the key piece of in- 
formation required to translate between egocentric and geocentric space is <fi. Such 
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translations will be required in any PI system which uses a combination of geocentric 
and egocentric vectors due to their different advantages, such as for HV updating 
versus homing. Without knowing <f), a point in egocentric space corresponds to a 
circle of possible geocentric locations, and a point in geocentric space corresponds 
to a circle of possible egocentric locations. 



Conclusions 



It is clear from the literature reviewed in the current work that motivations 
for using any particular coordinate system for modelling PI varied considerably, 
particularly in the arthropod literature. Under the novel static/dynamic vectorial 
classification scheme, it can be seen that a number of models should be considered as 
variants of the same class, simplifying analyses. For example, a number of PI models 
have been publi shed in both the arthropod and mammalian literature which used a 



GC or GS HV (IBernardet et all. 120081: iBurgess et all. 120 07; Ha ferlach et all 12 007 



Kim and Hallam. 20 001; iMcNauehton et all. 2006: 



Mittelstaedt and Mittelstaedt. 1973 



Vickerstaff and Di Paold . 120051 : IWittmann and Schweglerl . Il995l ). Indeed, one of the 
earliest PI models ever published, the Mittelstaedt model, was a member of the GS 
class and was used as the basis for the equational models presented in the paper. 
Being GC itself, it assumed a pair of sine and cosine input signals from the com 
pass. Although similar signals have been found in nature ( ILewis and Kristan Jr 



19981 : Miller et all I1991I ). most of the other GS-allothetic models used an array of 



multiple inputs, including some where the shape of the response funct ion of each 
input was not sinusoidal. For instance, IWittmann and Schweglerl ( 119951 ) assumed a 
single peaked ac t ivatio n function spread around a circular array of neurons, whereas 
Bernardet et al.l (120081 ) assumed a single active neuron indicating direction. Thus, 
for the special case of a perpendicularly oriented, two valued compass input (a GC 
representation), strict sinusoidal functions are required for geometrically exact PI, 
whereas for a larger number the exact shape of the function is less important (see 
Appendix D). Despite the range of possible neural network realizations of the GS 
class, the GS HV is fundamentally accumulating distances travelled along a set of 
geocentrically fixed directions. For the reasons discussed earlier, (see also Table 7) it 
possessed relatively few computational disadvantages, but a number of advantages. 
Amongst these is the fact it is the most suitable coordinate system for HV updating 
when using an allothetic compass. Given the significa nt increase in r obustness t o 



noise gained from using an allothetic compass for PI (ICheung et all 120071 . 120081 ). 



this property alone seems to be strong reason to expect some sort of GS HV to be 
used by biological PI systems where an allothetic compass is available. It is notable 
that in the mammalian literature, where an abundance of in vivo electrophysiolog- 
ical data exists from behaving animals, the few published PI models belong to the 
GS class. 

Which kind of HV is most suitable when only idiothetic directional cues are avail- 
able? If it is assumed that idiothetic cues are present as information about the ani- 
mal's rotation rate, then the EC HV update equations ( Merkle et al. . 2006bl ) appear 
to be simplest to implement (Table 3, EC HV Updating), followed by EP, although 
it is also possible to integrate rotation rate to estimate absolute heading for use 
with either geocentric PI systems (Table 4). However, the disastrous consequences 
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of PI using idiothetic directional cues are now well understood (ICheung et al.l . 120071 . 



20081 ). so extreme care must be taken during experimentation and analysis before 



concluding that any animal is navigating purely by PI in the complete absence of 
allothetic directional cues. 

Separating the HV into a directional and distance component seems to bring 
unnecessary complications, which apply to both egocentric and geocentric systems. 
The most obvious drawback from the equations for HV updating is the potential 
division-by-zero when the animal has a zeroed HV, but more generally the rate of 
change of the angular part of the HV increases towards infinity the closer the animal 
is to home. An EC HV has an increasing rate of change the further the animal is from 
home since the required response of the HV to rotation of the animal is proportional 
to the distance from home. Given the high rotation rates observed in dese rt ants, up 
to 4000° per second in the fastest species ( 1 Wehner and Srinivasanl . Il98ll ). a neural 
network implementation needs to have a neurophysiogically plausible method of 
coping with this phenomenon. In contrast, the GC HV has a maximum rate of 
HV change, limited simply by its maximum forward speed. Furthermore, GC HV 
coordinates can each be independently updated using information from only speed 
and directional cues, whereas for GP, EP and EC HVs, correct updating of each HV 
value requires information about speed, direction and also the other member of the 
pair of HV values. For these reasons a GC or GS HV representation seems to be a 
simpler and there fore more robust de s ign, w ith fewer potential sources of error. 

The model of iMiiller and Wehnerl (119881 ) employed a GP model to explain sys- 
tematic homing errors in desert ants, but failed to address the possibility that a 
similar model using a different coordinate system might produce the same errors. 
This paper has highlighted how, for the leaky integration model at least, identical 
results can be obtained regardless of the type of HV used. Whole classes of possible 
systematic error models may share this property, hence a more rigorous treatment 
of the subject appears necessary before attempting to use such errors as a probe of 
the internal workings of biological PI systems. 

While egocentric coordinates appear to be necessary and sufficient to control 
both homing and searching, it is impossible to exclude the possibility that biological 
PI systems may use a geocentric HV system, dynamically translating the HV into 
egocentric terms for controlling behaviour. This possibility was highlighted in the 
earlier description of the pendulum search model where analysis of the behaviour 
appeared consistent with the idea that the simulated animal was using a rotation 
rate sensor and an EC or EP HV during search, when in fact it used a GC system 
throughout. It demonstrated the difficulty of deducing the type of PI system being 
used based on external observations alone, but still demonstrated the utility of 
egocentric visualisation of homing and searching data. Thus it is prudent to consider 
using a variety of reference frames in studying PI behaviour, but care must be taken 
in drawing conclusions from the analyses. 

This paper has largely ignored the effects of noise and random errors on PI. 
Under this simplifying assumption, we have shown the extent of the mathematical 
equivalence of alternative coordinate systems for PI, and noted how difficult it may 
be to determine the internals of a PI system from external behaviour. We introduced 
a coordinate classification system which adequately covered all exact PI models 
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reviewed. Using this classification scheme, we showed that important properties 
of PI, especially related to HV updating, behave in similar ways for all members 
of each extended class. Therefore, this should form a useful foundation for future 
analyses of PI models, including the study of the effect of neural noise. 

Valid arguments were found in favour of and against using each coordinate sys- 
tem. On balance, the current theoretical analysis favours a geocentric and Cartesian 
or static vectorial HV when an allothetic compass cue is available. Future work is 
required to investigate the robustness of the different classes of HV in the face of 
random noise. This approach may offer the best hope for determining the type of PI 
coordinate system used in nature. We believe this will be an important step towards 
an understanding of the neurobiology of animal navigation. 

Appendix A: Derivation of HV Update Equations 

GP from GC 

Starting from the definition of GC HV updating (Table 3), and using the follow- 
ing relations (see Table 1): 



x = r cos 9 

y = r sin 9 (A.l) 

r = \/ x 2 + y 2 

it can easily be shown that r = s cos((f)—6) by differentiating r = ^ x 2 + y 2 . Defining 
z = x 2 + y 2 , such that r = yfz we find that z = 2xx + 2yy and r = using the 
chain rule. Using the relations given by Eqns IA.14 z expands to: 



z — 2 • r cos 9 ■ s cos (p + 2 ■ rsinO ■ s sin < 
z = 2rs(cos 9 cos <fi + sin 9 sin 0) 
z = 2rs cos(# — 0) 
z = 2rs cos(0 — 9) 

2rs cos((p — 9) 



2^~z 

r = — 2rs cos(0 — 9) 
2r KY ' 

r = s cos(0 — 9) 

Next we show that 9 = - sin(</> — 9) by differentiating x = r cos 9, giving x 
r cos 9 — r9 sin 9 by the product rule. Rewriting in terms of r9 sin 9 gives: 
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r cos 9 — x 

s cos(0 — 9) cos 9 — x 
s cos(0 — 9) cos 9 — s cos 
cos(0 — 9) cos 9 — cos 

(cos cos + sin sin 0) cos 9 — cos 

cos cos 2 + sin sin cos — cos 

cos 0(cos 2 9 — 1) + sin sin cos 

sin sin cos 9 — cos sin 2 9 

- (sin cos 9 — cos sin 9) 

r 

- sin(0 - 9) 

r 

EC from GC 

If = 0, the direction of the x and x' axes are aligned, and those of the y and y' 
axes are also aligned, giving the relation x' = —x and y' = —y. Rotation of the agent 
ACW by radians rotates the egocentric axes ACW relative to the geocentric axes, 
and therefore has the effect of rotating points in egocentric space CW by radians. 
Therefore (x',y') can be obtained by rotating (—x, —y) CW by radians, or ACW 
by —0 radians (from the standard matrix operation of rotation anticlockwise): 

x' = (—x) cos(— 0) — (— y) sin(— 0) 
y ' = (- x ) s in(-0) + (-y) cos(-0) 
x = —x cos — y sin 
y' = x sin — y cos 

Expressed in reverse (x,y) is obtained by rotating (—x 1 , —y') ACW by 0: 

x = (—x 1 ) cos(0) — (— y') sin(0) 
V = sin(0) + {-y') cos(0) 

x = —x' cos + y' sin 
y = —x 1 sin — y' cos 

Differentiating x', using the product rule, and substituting for x and y: 



r9 sin = 
r9 sin = 
r9 sin = 

-0sin0 = 

s 

-9sin9 = 

s 

-9sin9 = 

s 

-9sin9 = 

s 

-9sm9 = 

s 

9 = 
9 = 
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x = x<fi sin — x cos — y<p cos — y sin 

x = x<p sin — (s cos 0) cos <f) — y<fi cos — (s sin 0) sin 

i:' = x0 sin — ycj) cos — s(cos 2 + sin 2 0) 

i;' = x<p sin — y<p cos — s 

i/ = 0(x sin — y cos 0) — s 

x' = <py' — s 

Differentiating y f , using the product rule, and substituting for x and y: 



y' = x<fi cos + x sin + y0 sin — ?/ cos 

= x0 cos + (s cos 0) sin + y0 sin — (s sin 0) cos 
?/' = xcj) cos + y<p sin + s(cos sin — sin cos 0) 
y = x<p cos + ycj) sin 
y' = 0(x cos + y sin 0) 
2/' = -<K 

£P /rom GP 

It is clear that r' — r at all times, since both measure the distance between the 
animal and the home location. If = the x and x' axes are parallel and it is clear 
that 9' = 9 ± 7r. Rotating the animal by ACW will reduce by 0, but leave # 
unchanged. Therefore, the relationship 9' = 9 ± n — (or # = ± 7r + ) holds. 
Taking r = s cos(0 — 



Taking 0' = 9 ± vr - 0: 



s cos(0 — yj 

SCOs(0- (#'±7T + 0)) 

s cos(— ± 7r) 

S COs(#' ± 7r) 
— s COS 0' 
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e' = e - j> 

9' = - sin(0 -6)-<j> 
r 

9' = 4 sin(0 - 0) - 
r 

9' = ^sin(0- (0'±tt + 0)) -0 

0' = 4 sin(-0' ± tt) - 

r 

0' = -sin0'-0 

r 

Appendix B: Leaky integrator model 

Motion with a Fixed Velocity 



x 



dx 
~dt 



s cos0 — k D x .' 



dx 



XQ scos<f)/k D + x 
■. In (— s cos 0/ &d + x) \H = 
\x t = s cos 0/ &d + (#o — s cos 0/ e~ fel3< 



—k D / (it 
Jo 

-k D t 



d y ■ A\ b f yt d JL V \ M 
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Direct homing 



dr' , . r dr' , f* , 

— = -s - for .-. / = -fo / dt 
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.\\n(s/k D + r') |°, = —hot 
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Appendix C: Exact GS HV Updating Equations 

Any two non-parallel vectors, u, v, can form the basis of a GS coordinate system. 
Assuming they are both of unitary length, they can be defined by the direction the 
vectors point, 9 U ,9 V , giving vectors of the form: 
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u 



V = 



cos6 u 
sin 9 U 

cos 9 V 
sin 9 V 



A HV using this system, (u,v), indicates a GC position of: 



u 



cos 9 U 
sin 6L 



+ w 



cos ^ 
sin 6L 



To find the HV update equations for this system the following must be solved 
(from the GC HV Update equations given in Table 3): 



= u 



cos 9 U 
sin 9„ 



+ v 



cos 9 V 
sin 9 V 



= s 



cos <p 
sin0 



Assuming it to take the form a u cos 4> + b u sin and ?) to take the form a v cos + 
fe^sin^, it is straightforward to show that: 



sin^ 



a„ = 



sm(9 v - 6 U ) 

cos 6 V 
sm(9 v - 9 U ) 

sin 9 U 



sm(9 u - 9 V ) 

cos 9 U 
sin(6' u - 6>„) 



Giving the following result: 



sin 



u — s 



V = s 



s\n{9 v - 9 U ) 
sin(6' M - 4>) 
sin(9 u - 9 V ) 



Thus the update equation for each vector is a scaled and phase shifted sinewave. 
Rather than the maximum rate of change necessarily being in the direction of the 
coordinate's associated vector, as might have been expected, instead the rate of 
change is always zero in the direction of the other coordinate's vector (i.e. u — 
when = 9 V and vice versa). It is clear that the equations cannot be used if 
sin(^„ — 9 U ) = 0, confirming that only a pair of non-parallel static vectors can be 
used. 
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To construct a system using n vectors, each vector may form a pair with the 
previous vector, and a second pair with the following vector (including the last and 
first vector pairing), giving n pairs in total. For a set of vectors: 



cos 6i 
sin Qi 



forming the basis of an n valued HV (r±, • • • ,r n ), where the GC coordinate is 
found using: 



the following update equation may be used (where n + 1 = 1 and 1 — 1 = n): 



This gives a simple, geometrically exact HV updating system for any n > 2, 
provided only that no vectors adjacent in the list are parallel. 

Appendix D: Asymptotically Exact GS HV Update Equations 

For GS vectorial representations, a simple updating mechanism can be imple- 
mented using a transfer function / which, like a tuning curve, has a variable response 
which depends on the current heading <fi. The HV representation requires a ring of 
n neural units, each with the same transfer function /, but offset by — radians from 
the adjacent one. Thus fj(<f>) = fj+i((/> + ^f)- Then the HV update equation for 
axis j is fj = sfj(</>). 

Clearly, when n = 4 and fi = cos then f 2 = cos(0+|) = sin 0, / 3 = sin(0+|) = 
— cos0, and f\ = — cos(0 + f ) = — sin0, and is equivalent to the standard GC 
coordinate system, but with negative axes considered as distinct representations in 
their preferred direction. 

When n is small, the simple scheme presented gives an exact HV only for spe- 
cial cases of /, like the cosine/sine functions when n — 4. However, as n — > 
oo, the structure approaches a continuous representation in angular space hence 
lim 27r / n ^.o fo{4>) = fe{4> + 0)- Below it will be shown that the exact shape of function 
/ becomes unimportant for large n. 

For any unbiased function / e.g. symmetrical about 0, then re = sf $((/)). For an 
unbiased / we have: 






and 
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where w is a constant. 



/TV 
fo(9 + 0) sin(fl + <f))d9 = 
-7T 

/7T 
Xg COS 6* 6?$ 
-7T 

/7T 
cos 
-7T 

= s / f (<p + 9) cos 9d9 



From the definitions of an unbiased / we have: 



/7T 
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■IT 

f ((f) + 9) cos (f) cos 9 d9 — s / / o (0 + 0) sin0sin0d0 
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f o(4> + 9) cos 9 d9 + s cos (f) / f ((f) + 9) sin 9 d9 
IT J —IT 

/ o (0 + 0) cos 9d9 + ssin0 / / o (0 + 9) sin 9d9 

IT J —TV 



Summing the two results we get: 



/7T 
/o(0 + 9) cos 9d9 
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/7T 
/o(0 + 9) cos 9d9 
-IT 

/7T 
fo((f> + 9) cos 9d9 
-IT 



i(0) 
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Similarly, y((f>) = wsin<f). Therefore as n — > oo, the simple GS HV update 
equation r e = sfe((p) results in the same HV representation as the standard GC 
equational model. However, in the general GS case, / may be of any form as long 
as it is unbiased. 

These two methods of updating a GS HV prove that varying the number of static 
vectors cannot be assumed to yield different PI behaviours, since they can be made 
to reproduce exactly a GC HV. At least two sets of update equations exist (shown 
above) which make all noise- free static vectorial representations equivalent. 
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Figure 1: Four ways to represent the same spatial relationship between animal and home. 'A' is 
the animal's location, the associated arrow shows the orientation of its body axis. 'H' shows the 
home location. Shown are HVs for each of the four 'standard' coordinate systems considered in this 
paper, (a) geocentric Cartesian (GC), (b) geocentric polar (GP), (c) egocentric Cartesian (EC), 
(d) egocentric polar (EP). The angle <j> (shown only for the GC case, but relevant for all cases) is 
the absolute compass heading, and is not part of the HV. The angles 6, 9' and <p are measured 
positive in the anti-clockwise direction in radians. 
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Figure 2: The leaky integration model of systematic homing errors fitted to data from the dog- leg 
experiments of iMiiller and Wehnen (l988) where the first section was 10m and the second section 
5m in length, x-axis: angle turned by the ant between the first and second section (clockwise from 
direction of first section), y-axis: homing angle after completing the second section (clockwise from 
direction opposite that of first section). Solid line: predicted homing angle with a decay rate of 
ko = 0.0185s -1 and a walking speed of 0.33ms -1 . Dashed line: predicted homing angle with decay 
rate of zero, i.e. the geometrically correct homing angle. Stars: data points from the experiments. 
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Figure 3: The l eaky integration model of systematic homing errors fitted to data from the dog- leg 
experiments of iMiiller and Wehnerl (Il988h where the first section was 5m and the second section 
10m in length, x-axis: angle turned by the ant between the first and second section (clockwise from 
direction of first section), y-axis: homing angle after completing the second section (clockwise from 
direction opposite that of first section). Solid line: predicted homing angle with a decay rate of 
ko = 0.0185s -1 and a walking speed of 0.33ms -1 . Dashed line: predicted homing angle with decay 
rate of zero, i.e. the geometrically correct homing angle. Stars: data points from the experiments. 
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Figure 4: The leaky integration model of systematic homing errors fitted to data from the dog-leg 
experiments of iMiiller and Wehnen dl988h where the first section was 10m and the second section 
10m in length, x-axis: angle turned by the ant between the first and second section (clockwise from 
direction of first section), y-axis: homing angle after completing the second section (clockwise from 
direction opposite that of first section). Solid line: predicted homing angle with a decay rate of 
ko = 0.0185s -1 and a walking speed of 0.33ms -1 . Dashed line: predicted homing angle with decay 
rate of zero, i.e. the geometrically correct homing angle. Stars: data points from the experiments. 
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Figure 5: The leaky integratio n model of systematic homin g errors fitted to data from the long 
straight journey experiments of Sommer and Wehnerl t|2004h . x-axis: distance walked on outward 
leg. y-axis: distance walked on return leg. Solid line: predicted homing distance with a decay rate 
of kr> = 0. 00171s -1 and a walking speed of 0.33ms -1 . The geometrically correct homing distance 



is equal to the outward distance. Stars: data points from the experiments. 
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Figure 6: Homing trajectory generated by the searching model after an L shaped outwards excur- 
sion. The animal begins at location (0,0). The outward journey consists of one leg 10m to the 
right, followed by 5m upwards. Homing begins immediately after the second leg ends, and is gen- 
erated by the searching model (including the initial, straight homing leg and the subsequent search 
pattern centred on the home location). The animal's forward speed is fixed at 1ms -1 throughout 
the entire journey. The model parameters used (see Table 3, Searching) are s = 1, k\ = 2.7973 and 
k 2 = 1.308. The simulation covers 400 seconds. The trajectory was produced by numerical inte- 
gration of the HV update and searching equations using GC coordinates, using the Runge-Kutta 
4th order method with a fixed step size of 0.01 seconds. No noise is present in the simulation. 
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Figure 7: Search density profile of a single search pattern lasting for 10 6 seconds, using search 
parameters s = 1, k\ = 2.7973 and fc 2 = 1.308 (see Table 3, Searching). The animal was returning 
from an outward foraging journey as shown in Fig. [5] The trajectory was produced by numerical 
integration of the HV update and searching equations using GC coordinates, using the Runge-Kutta 
4th order method with a fixed step size of 0.01 seconds. Search density was recorded (excluding 
the outward part of the journey) on an 800 by 800 grid covering an area of 8 by 8 metres centred on 
the home location. The plot shows a slice through the two dimensional probability density at y=0. 
For comparison, the dotted line shows a fitted radially symmetrical bivariate normal distribution. 
The fitted value of a was 0.90028m assuming a correlation of p = (i.e. the fitted curve is for a 
radially symmetrical bivariate normal). 
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Figure 8: Plot of a single search pattern lasting for 400 seconds, showing the EC HV and rotation 
rate, using search parameters s = 1, ki = 2.7973 and ki — 1.308 (see Table 3, Searching). The 
animal was returning from an outward foraging journey as shown in Fig. [6l The trajectory was 
produced by numerical integration of the HV update and searching equations using EC coordinates, 
using the Runge-Kutta 4th order method with a fixed step size of 0.01 seconds. The two main 
sets of loops in the plot correspond to the parts of the geocentric search when the ant is circling 
towards the nest in an anticlockwise (positive rotation rate) or a clockwise (negative rotation rate) 
direction. 
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From \To 


GC 


GP 


EC 


EP 


GC 




r = \J x 2 + y 2 
9 = atan2(y, x) 


x = —x cos — y sin (f) 
y' = x sin 4> — y cos </> 


r — \/ x 2 + y 2 
9' = atan2(y, x) ± 7r — (f> 


GP 


x = r cos 9 
y = r sin 9 




x — — rcos(0 — 9) 
y' — r sin((/> — 9) 


r = r 
9' = 9 ± tt - 4> 


EC 


x = —x' cos (f> + y sin <j> 
y = —x' sin <j> — y' cos <j> 


r = \J x' 2 + y' 2 
9 = atan2(y',x') ± tt + <p 




r = \J x' 2 + y' 2 
9' = atan2(y',x') 


EP 


x = —r cos(9' + <p) 
y = —r' sin(#' + <p) 


r = r 

= 6' ± 7T + <)) 


x = r cos 9' 
y' = r' sin 9' 





Tabic 1: Equations for converting home vectors between four standard coordinate systems: geocentric Cartesian (GC), geocentric polar (GP), egocentric 
Cartesian (EC) and egocentric polar (EP). The equations convert from the system named in the left column to the system named in the top row. (f> 
indicates the animal's compass heading measured anti-clockwise from the direction of the x axis in radians. atan2(y,x) is defined as the angle from the x 
axis to a line from the origin passing through the point (x,y), again measured positive anti-clockwise, and can be calculated using arctan^ + kir, where 
fc = 0ifx>0, fc = lifx<0 and y > 0, and fc = — 1 if x < and y < 0. 



Usage 


Abbreviation 


Meaning 


General 


PI 


Path integration 


HV 


Home vector 


Standard coordinate systems 


GC 


Geocentric Cartesian 


GP 


Geocentric polar 


EC 


Egocentric Cartesian 


EP 


Egocentric polar 


Extended coordinate systems 
(vectorial basis representation) 


GS 


Geocentric static 


GD 


Geocentric dynamic 


ES 


Egocentric static 


ED 


Egocentric dynamic 


Path integration (PI) inputs 


s 


Speed (ms _1 ) 


4> 


Absolute heading (rad) 


Home vector (HV) values 




Geocentric x distance, egocentric x distance (rostrocaudal axis) (m) 


y, y' 


Geocentric y distance, egocentric y distance (left-right axis) (m) 


r, r' 


Geocentric radial HV distance, egocentric radial HV distance (m) 


8,8' 


Geocentric HV angle, egocentric HV angle (m) 


Constants / coefficients 


k s 


Speed coefficient during steering (s _1 ) 


k D 


Decay rate of first order decay of HV distance representations (leaky integrator) (s _1 ) 


k$ 


Coefficient of turning during steering (m _1 s _1 ) 


kn 


Search radius within which agent spends half of its time (m) 


k\,k2 


Weights egocentric y position of home and rotation rate 

to dynamically change homing/search path curvature (m — -'-s^^s -1 ) 



Table 2: Abbreviations used in this paper. SI units are included here for convenience but do not affect the form of the equations used in this work. 





HV Updat- 
ing 


Systematic 
Errors 


Homing 


Searching 


GC 


X — S COS (j) 

y = s sin 4> 


x = s cos <f> — Udx 
y = ssin</>- k D y 


4> — fc$ (x sin (j) — y cos <f>) 
s = ks(—x cos <p — y sin <j>) 


<j> = fa {x sin <j) — y cos <f>) — ki$ 


GP 


r = s cos(</> — 9) 

6 = - siaU - 9) 
r 


r = s cos(0 — 9) — k^r 

9=- sinU - 9) 
r 


4> = fc$r sin(</> — 9) 
s = —ksr cos(</> — 9) 


(f) = kirsm(<f> — 9) — k 2 <p 


EC 


x = <py' — s 

y = -<j>x' 


x = 4>y' — s — kpx' 
il = -4>x' - k D y 


4> = fc$y' 
s = ksx' 


4> = fay' - fa4> 


EP 


r' = — s cos 9 1 

9> = -sin 6' - 4> 

r 


r' = — s cos 9' — kur' 

9' = -t-sin^'-^ 

r 


4> = fc$r' sin#' 
s = ksr' cos 9' 


(j) = far' sa\9' — fc 2 



Table 3: Continuous-time equations in each of the standard coordinate systems describing HV updating, systematic PI errors, homing and searching. 
Dot notation indicates the first derivative with respect to time (for example x indicates two dots indicates the second derivative with respect to 
time. Symbols used: (x,y), (r,9), (x',y') and (r',9') are the GC, GP, EC and EP home vector respectively; s(ms _1 ), the animal's forward speed; </>, the 
animal's compass heading measured anticlockwise from the x axis direction in radians; fc/^s -1 ), a positive decay rate; /c$(s _1 m. _1 ), a positive turning 
rate coefficient; fcs^ 1 ), a positive forward speed coefficient; fci(s _2 m _1 ), ^(s -1 ), two positive constants controlling searching behaviour. Forward speed 
is held at a constant value during searching, s = 1, ki = 2.7973, k 2 = 1.308 gives a search pattern with a walking speed of I m s _1 and a search pattern 
which spends approximately half of the time within a radius of 1 m of the origin. To scale to arbitrary walking speed s and search radius fcij(m) use 
fa = 2.7973^-, fa = 1.3084. 





Allothetic directional cue 


Idiothetic directional cue 


GC 






Continuous 








x = s cos <j> 
y = s sin <j> 


X — S COS 1 / (pdt J 

/ r , \ 
y = s sin 1 / <pdt I 


GC 






Discrete 








%t+At —xt = S t + At COS (j) t+ At 


xt+At - xt = s t+ At cos (^2^ A0J 




Vt+At —yt = s t+ Atsin4> t+ At 


Vt+At -Vt = s t+A t sin 1 2 / A0I 


GP 






Continuous 








r = s cos ((ft — 6) 


r = s cos 1 1 / 0ai 1 — f 1 




= - sin (0-0) 

V 


s f f f ■ \ \ 

e= r sin {(j tdtj-e) 


GP 






Discrete 








n+At — n « st+At cos (4>t+At - St) 


rt+At - n ~ St + At cos A0 - e «) 




St+At - 0t ~ * +A * sin(t/) t+ Ai - 0t) 


0t +At fl*» St+At sin(2^A0 t ) 




n+At = y/ r t + S t+At + 2rts t+ AiCos (0t+Ai - ^i) 


»*i+At = \/ r t + S t+At + 2r * s t+At COS (A0 t+ At) 




ft+Ai - 0t = sin -1 ( St+At sin {</> t +At ~ <Pt) ) 
\rt+At J 


f t+Ai - 0< = sin~ 1 ( St+At s i n ( A</> t+ At) ) 
V»"t+Ai / 



Table 4: HV updating equations for GC and GP coordinate systems presented in continuous- and discrete-time, and for the case of allothetic or idiothetic 
directional cues. Symbol conventions are as for Table [3] Acj>t+At = <frt+At — 4>t- Using an absolute directional cue (true compass), rotation rate can only 
be estimate e.g. = lim At ^ . 



EC 

Continuous 



EC 

Discrete 



EP 

Continuous 



EP 

Discrete 



Allothetic directional cue 



c)l 



y - s 



at 



X 't+At - x 't ~ (0t+At - 4>t) y't - St+At 

y't+At -y't - - {4>t+At - 4>t)x' t 

4+At = x 't cos (0t+At - <£t) + y't sin (</>t+At - (fit) - st+At 
y't+At = ~ x 't sin (</>t+At - <t>t)+y' t COS ((f>t+At - 4>t) 



-s cos 8' 

at 



' t+At 



St+At cos ti t 



St+At ■ D l i i i \ 
— sm6» t - {(pt+At - <t>t) 



't+At 



y t+At 



( St+At . 

— sin {(pt+At 

\ r t+At 



S t+At + ^ r 't s t+At COS (<t>t+At - <t>t) 
) ) + 7T- {fa+At - 4>t) 



Idiothetic directional cue 



x' = <j>y' — s 

■ i ] i 

y = -<px 



X 't+At ~ x 't ~ ^t+Aty't ~ St+At 

y't+At - y't - -&<i>t+At x 't 

X 't+At = x 't cos (A<j>t+At) + y't sin (A<f> t +At) ~ s t +At 
y't+At = ~ x 't sin (A</>t+At) + y't cos (A(f> t +At) 



-s cos 9' 



r t+At r t 
„/ . S t +At 



~ -St+At cos tf t 
sm6' t - A(p t +At 



't+At 



\J r't + s^ +At + 2r' t s t +At cos (A<p t +At) 
sin (A0 t+At ) ) +tt- A<f> t +At 



-1 f St+At 



\ r' 



Table 5: HV updating equations for EC and EP coordinate systems presented in continuous and discrete time, and for the case of allothetic or idiothetic 
directional cues. Symbol conventions are as for Table [3l A</> t+ At = (frt+At — Using an absolute directional cue (true compass), rotation rate can only 
be estimate e.g. ^ = liniA 



At 





Homing 


Searching 


GC Discrete 


4>t+At -<t>t = k$ Ot sin 4> t - yt cos <j> t ) At 
st+At = —ks (xt cos <p t + yt sincA t ) At 


(pt+At 




= ki (x t sincAt - y t coscpt) - (f> t -At (l ~ + 2 <Pt 


GP Discrete 


<Pt+At — 4>t = k<s>n sin (4> t - 9 t ) At 
st+At = -ks'ft cos (4> t - t ) At 


<Pt+At 




= km sin (4> t - t ) - (pt-At (i - + 2( Pt 


EC Discrete 


4>t+At — 4>t = k§y' t At 
st+At = k s x' t At 


(pt+At 




= k lV 't-<Pt-At (l-^)+2& 


EP Discrete 


<t>t+At - 4>t = fc*r{ sin (e' t ) At 
st+At = k s r' t cos (e' t ) At 


(pt+At 




= kir't sine' t - <P t -At (l - ^) + 24>t 



Table 6: The generic homing equations and pendulum search model expressed as discrete time difference equations. Continuous to discrete approximation: 

J. ~ At AT ~ At At _ <Pt + At+<Pt-At-2<Pt 



Coordinate System 


GC 


CP 


EC 


EP 


Rate of HV change approaches infinity near home 




• 




• 


Rate of HV change increases with distance from home 






• 




HV updating needed in absence of physical translation 






• 


• 


Change in HV depends on current HV state 




• 


• 


• 


Further processing required to use allothetic directional cues 






• 


• 


Further processing required to use idiothetic directional cues 


• 


• 






HV alone insufficient to memorise and return to a geocentric location 






• 


• 


HV can be used for steering towards home 


o 


o 


o 


o 


Linear transform of HV sufficient for homing 






o 




HV can be used for searching for home (e.g. pendulum model) 


o 


o 


o 


o 


Exact update equations exist 


o 


o 


o 


o 


Can accomodate variable speeds of motion 


o 


o 


o 


o 


Can generate appropriate deviations from geometric PI (e.g. leaky integrator) 


o 


o 


o 


o 



Table 7: Comparison of the general properties of the main PI related equational models presented in this paper (see Table |3j). Filled circles (•), likely 
disadvantage; open circles (o), likely advantage. 



